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During  the  past  year,  we  have  significantly  furthered  our  studies  towards  meeting  the  goals  set  forth  in  the  funded 
proposal.  An  exciting  phenotype  has  been  identified  that  is  relevant  to  the  human  NF1  disease  and  to  the  potential 
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striking  learning  and  memory  defects  in  mutant  fish.  In  the  Look  lab,  the  zebrafish  models  of  the  NF1  tumor 
suppressor  linked  cancers  of  MPNST  and  glioma  have  been  completely  updated  to  make  them  optimal  for  the 
structure  function  studies.  Three  different  strategies  are  being  launched  to  further  accelerate  the  onset  of  these 
tumors  to  increase  the  precision  of  the  structure  function  studies  to  elucidate  the  mechanism  of  tumor  suppression 
by  NF1 .  We  are  making  arrangements  now  to  ship  our  NF1  zebrafish  lines  to  ZIRC,  so  they  are  available  to  our 
colleagues  in  the  NF1  and  zebrafish  communities.  These  genetically  engineered  lines  will  be  excellent  candidates 
for  use  in  screening  small  molecules  and  for  use  in  performing  genetic  screens  to  identify  modulators  of  this 
important  disease. 
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INTRODUCTION 

Neurofibromatosis  type  1  (NF1)  is  a  common  genetic  disorder  caused  by  mutations  in 
the  NF1  gene.  Patients  exhibit  pigmentation  abnormalities,  learning  disorders,  bone 
abnormalities  and  multiple  benign  and  malignant  tumors,  including  tumors  of  neural  crest 
origin  such  as  glioma  and  malignant  peripheral  nerve  sheath  tumors  (MPNSTs).  The 
NF1  gene  encodes  neurofibromin,  a  protein  of  over  4,000  amino  acids.  So  far,  only  a 
360  aa  region  of  neurofibromin  has  been  functionally  characterized  as  a  GAP-related 
domain  (GRD)  which  is  capable  of  accelerating  the  hydrolysis  of  GTP  bound  to  Ras, 
thus  down-regulating  Ras  activity.  Additional  important  functions  of  neurofibromin  are 
yet  to  be  discovered.  The  purpose  of  this  study  is  to  exploit  a  new  vertebrate  animal 
model  of  NF1.  We  capitalize  on  the  advantages  of  zebrafish  as  a  model  system  to 
address  pressing  questions  relevant  to  the  generation  of  new  and  effective  therapies  for 
NF1.  The  zebrafish  model  system  will  allow  the  Epstein  lab  in  Aim  1  to  perform  rapid  in 
vivo  rescue  “structure-function”  experiments  using  wild  type  and  mutant  NF1  genes,  as 
well  as  constructs  expressing  portions  of  the  NF1  protein,  in  order  to  determine  the 
functional  domains  of  neurofibromin  that  have  not  been  previously  identified  or  evaluated. 
These  structure-function  mutants  will  also  be  used  by  the  Look  lab  in  Aim  2  to 
understand  the  role  of  nfl  in  glioma  and  MPNST  tumorigenesis,  which  will  reveal 
clinically  relevant  mechanistic  information  about  neurofibromin  function  and  identify 
unrecognized  functional  domains  of  neurofibromin  appropriate  for  therapeutic  targeting. 


BODY 

Aim  1:  To  fully  characterize  nfla/b  compound  null  zebrafish  and  to  perform  structure- 
function  studies  of  Nfl  in  vivo  by  performing  rescue  experiments  in  znla/b  loss-of 
function  zebrafish  to  determine  which  aspects  of  the  mutant  phenotype  are  regulated 
by  GAP  and  non-GAP  domains  (Epstein  lab). 

The  Epstein  laboratory  has  continued  to  characterize  the  nfla/b  compound  null  zebrafish 
over  the  past  year.  An  exciting  phenotype  has  been  identified  that  is  relevant  to  the 
human  disease  and  to  the  potential  development  of  meaningful  therapies.  In 
collaboration  with  the  Granato  lab  at  Penn,  these  studies  have  identified  learning  and 
memory  defects  in  mutant  fish.  The  Granato  lab  has  developed  learning  and  memory 
assays  for  fish  ( 1 ,  2).  Briefly,  the  assays  utilize  the  fact  that  larval  fish  will  respond  to  a 
loud  noise  or  sudden  darkness  with  characteristic  body  movements.  These  movements 
can  be  recorded  with  a  high  resolution  video  camera  for  analysis.  After  repeated 
stimulations,  wild  type  fish  will  learn,  and  will  respond  with  greater  latency,  or  will  not 
respond  at  all.  After  an  hour  or  more  of  an  interruption  in  the  delivery  of  stimuli,  wild  type 
fish  will  remember,  and  will  respond  with  greater  latency  than  untrained  fish. 
Interestingly,  nfla/b  mutants  exhibit  both  learning  and  memory  defects  (Figure  1  and  2). 
Strikingly,  these  defects  can  be  corrected  by  the  addition  of  certain  small  molecules  to 
the  water.  Learning  can  be  corrected  by  drugs  (in  the  nanomolar  range)  that  modulate 
cAMP  (Figure  1),  while  memory  is  corrected  by  drugs  that  lower  active  Ras  signaling 
components,  including  PI3K  and  MAPK  inhibitors  (Figure  2).  Although  cAMP  signaling 
has  been  implicated  downstream  of  NF1  in  flies,  a  role  for  cAMP  in  NF1  pathology  in 
vertebrates  has  not  been  demonstrated.  Since  human  NF1  patients  have  learning  and 
memory  defects,  this  finding  may  have  direct  clinical  implications.  This  study  is  being 
prepared  for  publication. 
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NF1  mutant  larvae  show  deficit  in 
acoustic  startle  habituation 


Rolipram  rescues  NF1  startle  learning  deficit 


8-Br-cAMP  rescues  NF1  startle  learning  deficit 


Wortmannin/U0126  do  NOT  rescue  learning 


D=DMSO 

W=wortmannin  (1  uM| 
U=U0126  (1  uM) 


Figure  1:  nfla/b  mutant  fish  have  learning  defects  that  can  be  partially  corrected  by  Rolipram  and  8- 
Br-cAMP,  but  not  by  Wortmannin  or  U0126. 
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Wortmannin  rescues  Nfl  visual  memory  deficits 


8-Br-cAMP  does  NOT  rescue  memory  deficits 


500 


Memory  (Ras) 


Figure  2:  nfla/b  mutant  fish  have  memory  defects  that  can  be  partially  corrected  by  Wortmannin 
or  U0126,  but  not  by  8-Br-cAMP. 

The  Epstein  lab  has  also  made  progress  with  preparation  of  the  necessary  clones  for  the 
genetic  rescue  experiments  described  in  the  original  proposal.  Once  the  required 
plasmids  are  prepared,  we  will  be  ready  to  proceed  with  these  studies. 


Aim  2:  To  analyze  the  contributions  of  Nfl  GAP  and  other  functional  domains,  based  on 
Nfl  structure  function  analysis  in  Aim  1,  to  the  suppression  of  malignant  glioblastomas 
(CNS)  and  MPNSTs  (PNS),  which  we  have  shown  to  develop  in  conjunction  with  p53 
loss  (Look  lab). 

The  Look  lab  has  crossed  the  nfla/b  mutant  lines  with  the  p53  mutant  zebrafish  to 
generate  nfl a+/-;nf1  b-/-;p53-/-  fish  and  monitored  for  tumorigenesis.  The  current  nf1a+/- 
;nf1  b-/-;p53-/-  fish  line  began  to  develop  tumors  from  the  age  of  17  weeks  (Figure  3A). 
Before  36  weeks,  38%  nfl a+/-;nf1  b-/-;p53-/-  fish  developed  tumors,  including  melanoma 
(1/66,  indicated  by  the  black  arrow  in  Figure  3A),  high  grade  glioma  (3/66,  indicated  by 
the  red  arrows  in  Figure  3A;  an  example  fish  is  shown  in  Figure  3B)  and  MPNST  (21/66; 
an  example  fish  is  shown  in  Figure  3C),  whereas  no  tumor  was  observed  in  nf1b-/-;p53- 
/-  fish  (Figure  3A).  Within  one  year,  almost  all  nfl a+/-;nf1  b-/-;p53-/-  fish  developed 
tumors,  the  majority  of  which  are  MPNSTs,  proving  that  our  nf1/p53  mutant  zebrafish 
serve  as  a  robust  model  for  MPNST  studies.  To  improve  monitoring  tumor  initiation  and 
progression  in  vivo  in  the  nf1/p53  mutant  zebrafish  lines,  we  have  crossed  the  nf1/p53 
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mutant  fish  to  a  sox10:EGFP  line  (early  neural  crest  cells  are  labeled  with  GFP),  as 
soxlO  is  consistently  expressed  in  schwannian  tumors  (3,  4). 


nf1b-/-;p53-/-  (n=50) 

—  nf1a+/-;nf1b-/-;p53-/-(n=66) 


Figure  3:  (A)  Tumor  surveillance  results  for  nf1b-/-;p53-/-  (n=50)  and  nfla+/-;nf1b-/-;p53-/- 
(n=66)  zebrafish.  Tumor-free  survival  is  defined  as  the  time  in  weeks  from  date  of  birth  to  the  date 
of  tumor  onset.  25  nf1a+/-;nf1b-/-;p53-/-  fish  developed  tumors  by  the  age  of  36  weeks,  while  no 
nf1b-/-;p53-/-  fish  developed  tumors  yet.  The  black  arrow  indicates  one  melanoma  identified  at 
the  age  of  17  weeks.  The  red  arrows  indicate  three  high  grade  gliomas  identified  at  the  age  of  19 
and  27  weeks.  (B)  27-week-old  nf1a+/-;nf1b-/-;p53-/-  fish  with  high-grade  glioma  (arrow).  (C)  36- 
week-old  nfl  a+/-;nf1  b-/-;p53-/-  fish  with  MPNST  (arrow)  and  a  nevi  in  the  tail  fin  (arrowhead). 


Currently,  the  earliest  gliomas  developed  in  the  nfl a+/-;nf1  b-/-;p53-/-  line  were  observed 
from  19  weeks,  and  the  earliest  MPNSTs  were  observed  after  20  weeks.  To  improve 
the  model  for  the  nfl  structure-function  analysis  using  the  NF1  constructs  from  Aim  1, 
Shuning  He  in  the  Look  lab  is  taking  three  independent  strategies  to  further  accelerate 
the  tumor  onset  and  increase  tumor  penetrance. 

i.  Introducing  loss  of  pten  into  the  nf1/p53  mutant  lines.  Pten  (phosphatase  and  tensin 
homolog)  is  one  of  the  most  frequently  mutated  or  deleted  tumor  suppressors  in  human 
cancer.  Mutations  and  homozygous  deletions  of  pten  were  found  in  36%  of 
glioblastomas  (GBM),  the  most  malignant  subtype  of  glioma  (5).  Loss  or  decrease  of 
pten  expression  was  also  detected  in  a  majority  of  human  NFI-associated  MPNST 
lesions  (6).  There  are  two  orthologues  of  human  pten  in  the  zebrafish  genome,  namely 
ptena  and  ptenb  (7).  We  obtained  the  loss-of-function  pten  mutant  from  Dr.  Jeroen  den 
Hertog  and  bred  them  together  to  generate  a  ptena+/-;ptenb-/-  line.  Subsequently,  we 
established  a  nf1a+/-;nf1b+/-;p53+/-;ptena+/-;ptenb+/-  line  by  crossing  the  ptena+/- 
; ptenb-/-  fish  with  our  nfl a+/-;nf1  b-/-;p53-/-  fish.  We  have  incrossed  the  nf1a+/-;nf1b+/- 
;p53+/-;ptena+/-;ptenb+/-  line  to  obtain  all  possible  combinational  loss  of  nfl,  p53  and 
pten  in  our  zebrafish  model.  We  are  currently  performing  a  tumor  watch  on  these  fish  to 
determine  the  correlation  between  the  loss  of  tumor  suppressors  and  formation  of 
tumors  associated  with  neurofibromatosis  type  1.  Then  we  will  select  the  genotypes 
which  give  the  earliest  tumor  onset  and  high  tumor  penetrance  to  perform  the  structure- 
function  studies  using  the  NF1  constructs  from  Aim  1. 

ii.  Introducing  human  oncogenes  associated  with  GBM/MPNST  into  the  nfl a+/-;nf1  b-/- 

;p53-/-  mutant  line.  The  oncogenes  that  we  will  introduce  as  transgenes  are  mutant 
PDGFR  (8)  and  mutant  Met  (9),  which  are  frequently  mutated  in  human  GBM  and 
MPNST  (5,  8,  10,  11).  We  have  obtained  constructs  containing  constitutively  active 
mutants  of  human  PDGFR  from  Dr.  Eric  Holland  and  Met  from  Dr.  George  van  de 
Woude.  These  mutants  will  be  sub-cloned  downstream  of  the  zebrafish  soxlO  promoter. 
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Sox10:PDGFR  or  sox10:Met  will  be  injected  into  the  progeny  of  the  nfl  a+/-;nf1  b-/-;p53-/- 
\sox10:EGFP  zebrafish.  We  expect  that  injection  of  these  oncogenes  will  greatly 
accelerate  the  onset  of  high  grade  glioma  and  MPNST.  We  will  perform  a  weekly  tumor 
watch  to  determine  the  time  of  tumor  onset  in  the  injected  animal,  and  derive  stable 
transgenic  lines  for  the  nfl  structure-function  analysis.  We  will  select  the  construct  with 
the  best  GFP-expressing  glioma  acceleration  for  the  NF1  structure-function  studies, 
iii.  Introducing  mutations  in  H3F3A  and  ATRX  into  the  nf1a+/-;nf1b-/-;p53-/-  mutant  line. 
It  is  recently  found  that  H3F3A,  ATRX,  p53  and  nfl  are  the  most  frequently  mutated 
genes  in  pediatric  GBM  (72).  H3F3A  encodes  the  replication-independent  histone  3 
variant  H3.3  and  ATRX  (a-thalassaemia/mental  retardation  syndrome  X-linked)  is  a 
chromatin  remodeling  factor  required  for  H3.3  incorporation  at  pericentric 
heterochromatin  and  telomeres  (72).  The  H3.3/ATRX  perturbation  was  suggested  to 
play  a  central  role  in  pediatric  GBM.  The  Look  lab  is  applying  an  emerging  genome 
editing  technology  called  CRISPR-Cas  (73,  14)  to  incorporate  the  deficiencies  of  the 
H3F3A  and  ATRX  genes,  as  found  in  human  patients,  into  the  nfl a+/-;nf1  b-/-;p53-/- 
zebrafish.  We  expect  that  the  incorporation  of  H3F3A  and  ATRX  mutations  will  enhance 
high  grade  glioma  development  in  our  zebrafish  lines.  Once  established,  we  will  use  the 
NF1  constructs  from  Aim  1  to  perform  nfl  structure-function  analysis  in  these  zebrafish 
to  dissect  the  role  of  NF1  in  pediatric  high  grade  gliomas  associated  with 
neurofibromatosis  type  1. 


KEY  RESEARCH  ACCOMPLISHMENTS 

•  Characterization  of  learning  defects  in  nfl  mutant  zebrafish 

•  Characterization  of  memory  defects  in  nfl  mutant  zebrafish 

•  Preparation  of  NF1  constructs  for  the  structure-function  analysis 

•  Generation  of  nf1/p53/pten  mutant  zebrafish 

•  Preparation  of  glioma/MPNST  models  for  the  structure-function  analysis 


REPORTABLE  OUTCOMES 

•  Manuscript:  Zebrafish  neurofibromatosis  type  1  genes  have  redundant  functions 
in  tumorigenesis  and  embryonic  development.  Shin  J,  Padmanabhan  A,  de  Groh 
ED,  Lee  JS,  Haidar  S,  Dahlberg  S,  Guo  F,  He  S,  Wolman  MA,  Granato  M, 
Lawson  ND,  Wolfe  SA,  Kim  SH,  Solnica-Krezel  L,  Kanki  JP,  Ligon  KL,  Epstein 
JA,  Look  AT.  Disease  Model  and  Mechanism.  2012  Nov;5(6):881-94. 

•  Manuscript:  Evaluation  and  application  of  modularly  assembled  zinc-finger 
nucleases  in  zebrafish.  Zhu  C,  Smith  T,  McNulty  J,  Rayla  AL,  Lakshmanan  A, 
Siekmann  AF,  Buffardi  M,  Meng  X,  Shin  J,  Padmanabhan  A,  Cifuentes  D, 
Giraldez  AJ,  Look  AT,  Epstein  JA,  Lawson  ND,  Wolfe  SA.  Development  138(20): 
4555-64,2011. 

•  Abstract:  Zebrafish  neurofibromatosis  type  1  genes  have  redundant  functions  in 
tumorigenesis  and  embryonic  development.  Look  AT.  1st  Cold  Spring  Harbor 
Asia  conference  on  Fishing  for  Answers:  Zebrafish  Models  of  Human 
Development  and  Disease.  Suzhou,  China,  2012. 

•  Presentation:  Zebrafish  neurofibromatosis  type  1  genes  have  redundant 
functions  in  tumorigenesis  and  embryonic  development.  Look  AT.  1st  Cold 
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Spring  Harbor  Asia  conference  on  Fishing  for  Answers:  Zebrafish  Models  of 
Human  Development  and  Disease.  Suzhou,  China,  2012. 

•  Presentation:  Zebrafish  neurofibromatosis  type  1  genes  have  redundant 
functions  in  tumorigenesis  and  embryonic  development.  Look  AT.  SickKids 
Research  Institute,  Toronto,  Canada,  2012 

•  Employment  or  research  opportunities  applied  for  and/or  received  based  on 
experience/training  supported  by  this  award:  Jimann  Shin  (Ph.D.)  entered  a 
post-doc  training  program  in  Lila  Solnica-Krezel  Lab,  Washington  University  in 
St.  Louis. 

•  Employment  or  research  opportunities  applied  for  and/or  received  based  on 
experience/training  supported  by  this  award:  Arun  Padmanabhan  received 
MD/PhD  degree  from  the  Univ.  of  Pennsylvania,  2013 


CONCLUSION 

During  the  past  year,  we  have  significantly  furthered  our  studies  towards  meeting  the 
goals  set  forth  in  the  funded  proposal.  An  exciting  phenotype  has  been  identified  that  is 
relevant  to  the  human  NF1  disease  and  to  the  potential  development  of  meaningful 
therapies.  In  collaboration  with  the  Granato  lab  at  Penn,  the  Epstein  lab  has  discovered 
striking  learning  and  memory  defects  in  mutant  fish.  In  the  Look  lab,  the  zebrafish 
models  of  the  NF1  tumor  suppressor  linked  cancers  of  MPNST  and  glioma  have  been 
completely  updated  to  make  them  optimal  for  the  structure  function  studies.  Three 
different  strategies  are  being  launched  to  further  accelerate  the  onset  of  these  tumors  to 
increase  the  precision  of  the  structure  function  studies  to  elucidate  the  mechanism  of 
tumor  suppression  by  NF1.  We  are  making  arrangements  now  to  ship  our  NF1 
zebrafish  lines  to  ZIRC,  so  they  are  available  to  our  colleagues  in  the  NF1  and  zebrafish 
communities.  These  genetically  engineered  lines  will  be  excellent  candidates  for  use  in 
screening  small  molecules  and  for  use  in  performing  genetic  screens  to  identify 
modulators  of  this  important  disease. 
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Abstract:  Zebrafish  neurofibromatosis  type  1  genes  have  redundant  functions  in 
tumorigenesis  and  embryonic  development. 

Thomas  Look1,  Jimann  Shin1,  Arun  Padmanabhan2,  Eric  D  de  Groh2,  Jeong-Soo  Lee1, 
Jonathan  A  Epstein2 

^ana  Farber  Cancer  Institute,  Ped/Onc,  Boston,  MA,  02215 

2University  of  Pennsylvania,  Department  of  Cell  and  Developmental  Biology,  421  Curie 
Blvd,  PA,  19104 

Neurofibromatosis  type  1  (NF1)  is  a  common,  dominantly  inherited  genetic  disorder  that 
results  from  mutations  in  the  neurofibromin  1  (NF1)  gene.  Affected  individuals 
demonstrate  abnormalities  in  neural  crest-derived  tissues  that  include  hyperpigmented 
skin  lesions  and  benign  peripheral  nerve  sheath  tumors.  NF1  patients  also  have  a 
predisposition  to  malignancies  including  juvenile  myelomonocytic  leukemia  (JMML), 
optic  glioma,  glioblastoma,  schwannoma,  and  malignant  peripheral  nerve  sheath  tumors 
(MPNSTs).  In  an  effort  to  better  define  the  molecular  and  cellular  determinants  of  NF1 
disease  pathogenesis  in  vivo,  we  employed  targeted  mutagenesis  strategies  to  generate 
zebrafish  harboring  stable  germline  mutations  in  nfla  and  nflb,  orthologues  of  NF1. 
Animals  homozygous  for  loss-of-function  alleles  of  nfla  or  nflb  alone  are  phenotypically 
normal  and  viable.  Homozygous  loss  of  both  alleles  in  combination  generates  larval 
phenotypes  that  resemble  aspects  of  the  human  disease  and  results  in  larval  lethality 
between  7  and  10  days  post  fertilization,  nfl -null  larvae  demonstrate  significant  central 
and  peripheral  nervous  system  defects.  These  include  aberrant  proliferation  and 
differentiation  of  oligodendrocyte  progenitor  cells  (OPCs),  dysmorphic  myelin  sheaths, 
and  hyperplasia  of  Schwann  cells  and  sympathetic  neurons.  Loss  of  nfl  contributes  to 
tumorigenesis  as  demonstrated  by  an  accelerated  onset  and  increased  penetrance  of 
high-grade  gliomas  and  MPNSTs  in  adult  nf1a+/-;  nflb-/-;  p53e7/e7  animals.  Importantly, 
we  identify  and  quantitatively  analyze  a  novel  melanophore  phenotype  in  nfl -null  larvae, 
providing  the  first  animal  model  of  the  pathognomonic  pigmentation  lesions  of  NF1. 
Together,  these  findings  support  a  role  for  nfla  and  nflb  as  potent  tumor  suppressor 
genes  that  also  function  in  the  development  of  both  central  and  peripheral  glial  cells  as 
well  as  melanophores  in  zebrafish. 
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Zebrafish  neurofibromatosis  type  1  genes  have 
redundant  functions  in  tumorigenesis  and  embryonic 
development 

Jimann  Shin1'*'*,  Arun  Padmanabhan2'*,  Eric  D.  de  Groh2'*,  Jeong-Soo  Lee1§,  Sam  Haidar3,  Suzanne  Dahlberg1, 

Feng  Guo1,  Shuning  He1,  Marc  A.  Wolman2,  Michael  Granato2,  Nathan  D.  Lawson4,  Scot  A.  Wolfe4,  Seok-Hyung  Kim5, 

Lilianna  Solnica-Krezel6,  John  P.  Kanki1,  Keith  L.  Ligon3,  Jonathan  A.  Epstein2,11  and  A.  Thomas  Look1,11 

SUMMARY 

Neurofibromatosis  type  1  (NF1 )  is  a  common,  dominantly  inherited  genetic  disorder  that  results  from  mutations  in  the  neurofibromin  7  [NF1)  gene. 
Affected  individuals  demonstrate  abnormalities  in  neural-crest-derived  tissues  that  include  hyperpigmented  skin  lesions  and  benign  peripheral 
nerve  sheath  tumors.  NF1  patients  also  have  a  predisposition  to  malignancies  including  juvenile  myelomonocytic  leukemia  (JMML),  optic  glioma, 
glioblastoma,  schwannoma  and  malignant  peripheral  nerve  sheath  tumors  (MPNSTs).  In  an  effort  to  better  define  the  molecular  and  cellular 
determinants  of  NF1  disease  pathogenesis  in  vivo,  we  employed  targeted  mutagenesis  strategies  to  generate  zebrafish  harboring  stable  germline 
mutations  in  nfla  and  nflb,  orthologues  of  NFL  Animals  homozygous  for  loss-of-function  alleles  of  nfla  or  nflb  alone  are  phenotypically  normal 
and  viable.  Homozygous  loss  of  both  alleles  in  combination  generates  larval  phenotypes  that  resemble  aspects  of  the  human  disease  and  results 
in  larval  lethality  between  7  and  10  days  post  fertilization.  nf7-null  larvae  demonstrate  significant  central  and  peripheral  nervous  system  defects. 
These  include  aberrant  proliferation  and  differentiation  of  oligodendrocyte  progenitor  cells  (OPCs),  dysmorphic  myelin  sheaths  and  hyperplasia  of 
Schwann  cells.  Loss  of  nfl  contributes  to  tumorigenesis  as  demonstrated  by  an  accelerated  onset  and  increased  penetrance  of  high-grade  gliomas 
and  MPNSTs  in  adult  nf1a+/~;  nf1b~7~;  p53e7/e7  animals.  n/J-null  larvae  also  demonstrate  significant  motor  and  learning  defects.  Importantly,  we 
identify  and  quantitatively  analyze  a  novel  melanophore  phenotype  in  nf7-null  larvae,  providing  the  first  animal  model  of  the  pathognomonic 
pigmentation  lesions  of  NF1.  Together,  these  findings  support  a  role  for  nfla  and  nflb  as  potent  tumor  suppressor  genes  that  also  function  in  the 
development  of  both  central  and  peripheral  glial  cells  as  well  as  melanophores  in  zebrafish. 


INTRODUCTION 

Type  1  neurofibromatosis  (NF1)  is  an  autosomal  dominant 
inherited  genetic  disorder  characterized  by  pigmented  birthmarks 
known  as  cafe-au-lait  spots,  cutaneous  and  plexiform 
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neurofibromas  arising  in  the  glial  cells  of  the  peripheral  nervous 
system  (PNS),  optic  pathway  gliomas,  cardiovascular  abnormalities 
and  learning  defects  (Williams  et  al.,  2009).  The  disease  results  from 
mutations  in  the  NF1  gene,  encoding  the  large  protein 
neurofibromin,  which  contains  a  GTPase-activating  protein-related 
domain  (GRD)  capable  of  inactivating  the  RAS  proto-oncogene 
(Cawthon  et  al.,  1990;  Viskochil  et  al.,  1990;  Wallace  et  al.,  1990). 
Thus,  NF1  loss  results  in  aberrant  activation  of  Ras  signaling,  which 
may  predispose  NF1  patients  to  a  variety  of  cancers  (Cichowski 
and  Jacks,  2001).  Heterozygous  Nfl  mutant  mice  develop 
pheochromocytoma  and  myeloid  leukemia,  whereas  the 
conditional  loss  of  Nfl  in  a  p53- deficient  background  results  in 
highly  penetrant  malignant  astrocytoma  formation  (Jacks  et  al., 
1994;  Zhu  et  al.,  2005a;  Powers  et  al.,  2007).  Furthermore,  two  recent 
reports  have  identified  NF1  mutations  in  approximately  15-23%  of 
human  glioblastoma  patients  (Parsons  et  al.,  2008;  The  Cancer 
Genome  Atlas  Research  Network,  2008).  Although  these  studies 
demonstrate  a  strong  link  between  NF1  function  and  high-grade 
glioma,  the  crucial  signaling  pathways  governing  the  development 
of  tumorigenesis  remain  to  be  elucidated.  An  animal  model 
facilitating  the  rapid  interrogation  of  epistatic  and  functional 
relationships  within  signaling  pathways  would  serve  as  a  valuable 
tool  for  probing  the  pathology  underlying  NFl-induced  cell 
transformation. 

We  recently  developed  a  zebrafish  model  of  NF1  deficiency  using 
antisense  morpholino  oligonucleotides  to  produce  transient  gene 
knockdown  (Padmanabhan  et  al.,  2009;  Lee  et  al.,  2010).  Two 
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zebrafish  orthologues  were  identified  that  are  highly  homologous 
to  human  NF1  at  the  amino  acid  level,  sharing  approximately  84% 
identity,  including  91-93%  identity  within  the  GRD.  Both  genes 
maintain  syntenic  relationships  with  human  NF1  on  chromosome 
17qll.2  and  are  probably  the  result  of  the  well-described  genomic 
duplication  event  that  occurred  early  in  the  evolution  of  teleosts 
(Amores  et  al.,  1998).  In  our  previous  work  with  nfl  morphants, 
we  observed  defects  in  both  cardiovascular  and  nervous  system 
development.  However,  due  to  the  transient  nature  of  morpholino 
gene  knockdown,  the  analysis  of  nfl -deficient  phenotypes  beyond 
the  first  3  days  of  life  was  not  possible. 

We  report  here  the  generation  of  stable  mutant  nfl  zebrafish 
lines,  using  both  zinc  finger  nuclease  (ZFN)  and  targeting  induced 
local  lesions  in  genomes  (TILLING)  strategies,  and  the  detailed 
phenotypic  analysis  of  this  new  animal  model  of  human  NF1.  We 
have  successfully  generated  several  independent  null  alleles  of  nfla 
and  nfl b.  Mutant  larvae  carrying  at  least  one  wild-type  nfla  or 
nflb  allele  are  viable,  fertile  and  show  no  obvious  phenotypes  during 
early  development.  By  contrast,  nfla”1”;  nflb  larvae  exhibit  overt 
pigmentation  defects  as  early  as  6  days  post  fertilization  (dpf )  and 
do  not  survive  beyond  10  dpf.  Beginning  at  4  dpf,  nfla nflb ”r” 
larvae  exhibit  hyperplasia  of  oligodendrocyte  progenitor  cells 
(OPCs)  and  Schwann  cells,  as  well  as  melanophore  hypoplasia. 
Defects  resulting  from  the  loss  of  nfl  in  pigment  cell  and  glial  cell 
lineages  mirror  those  often  observed  in  the  tissues  of  NF1  patients. 
In  a p53  mutant  background  (p53e7le7),  nfla*1”;  nflb ”'”  fish  develop 
high-grade  gliomas  and  malignant  peripheral  nerve  sheath  tumors 
(MPNSTs),  demonstrating  a  tumor- suppressor  function  for  the 
zebrafish  nfl  orthologues.  Therefore,  we  have  developed  and 
characterized  a  heritable  zebrafish  model  of  NF1  that  exhibits 
clinical  hallmarks  of  the  disorder,  including  nervous  system  defects 
and  increased  susceptibility  to  tumorigenesis.  Furthermore,  nfl 
mutant  zebrafish  represent  the  first  vertebrate  model  of  the 
pathognomonic  pigmentation  lesions  associated  with  NF1. 

RESULTS 

Generation  of  zebrafish  nfla  and  nflb  mutants 

We  previously  identified  two  zebrafish  orthologues  of  human  NF1, 
nfla  and  nflb,  and  described  the  phenotypes  that  result  from  their 
loss  of  function  during  early  development  induced  by  antisense 
morpholino  oligonucleotides  (Padmanabhan  et  al.,  2009;  Lee  et  al., 
2010).  Although  this  technology  readily  permits  transient 
knockdown  of  gene  expression,  its  efficacy  is  limited  to  only  the 
first  few  days  of  life.  In  an  effort  to  gain  a  better  understanding  of 
the  roles  of  nfla  and  nflb  during  development,  as  well  as  in  cancer 
predisposition,  we  employed  multiple  approaches  to  develop  stable 
lines  of  zebrafish  harboring  germline  mutations  in  each  of  these 
genes.  Using  a  modular  approach  (Zhu  et  al.,  2011),  zinc  finger 
nucleases  (ZFNs)  were  engineered  with  binding  specificities  directed 
to  exon  26  of  nfla  and  exon  17  of  nflb  (Fig.  1A,B).  Paired  ZFN 
mRNAs  were  injected  into  zebrafish  embryos  and  independent 
target- specific  mutant  alleles  for  nfla  ( nfla A5  and  m/I«A8)  and  nflb 
(nflb*10  and  nflb A5S)  were  identified  in  the  FI  generation  (Fig.  1C; 
supplementary  material  Fig.  S1A-D).  Each  of  these  mutations 
included  a  deletion  and/or  insertion  within  a  coding  exon  that 
resulted  in  a  frameshift,  introducing  premature  stop  codons  that 
would  be  expected  to  truncate  the  neurofibromin  protein  upstream 
of  the  GRD  (Fig.  1D,E).  In  a  separate  effort,  we  screened  a  library 
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of  N-ethy\-N- nitrosourea  (ENU)-mutagenized  zebrafish  by  targeting 
induced  local  lesions  in  genomes  (TILLING)  (Wienholds  et  al.,  2002) 
and  identified  a  single  founder  harboring  a  nonsense  mutation  in 
exon  29  of  nfla  ( nflaL1247X )  (Fig.  1C;  supplementary  material  Fig. 
S1E-G).  To  confirm  that  the  targeted  alleles  disrupted  production 
of  full-length  protein,  we  performed  western  blots  using  an  antibody 
that  should  recognize  both  nfla  and  nflb  with  extracts  prepared 
from  3  dpf  wild-type,  «/l«A5/AS;  nflb*1*,  nfla*1*;  nflb*101*10  and 
«/lflA5/A5;  nflb*101*10  larvae  (Fig.  IF).  We  observed  a  complete  loss 
of  Nfl  signal  in  the  double-homozygous  null  extracts.  We  detected 
only  low  levels  of  protein  expression  in  w/2aA5/A5;  nflb*1*  mutant 
extracts  as  compared  with  wild-type  or  nfla*1*;  nflb*101*10  mutant 
extracts,  which  might  reflect  differences  in  expression  levels  of  the 
two  orthologues  at  3  dpf.  However,  we  cannot  rule  out  the  possibility 
that  the  neurofibromin  antibody  we  used  recognizes  the  two 
proteins  with  different  affinities.  We  generated  separate  zebrafish 
lines  with  distinct  null  alleles  of  both  nfla  and  nflb  to  provide 
evidence  that  the  observed  phenotypes  were  in  fact  due  to  nfl  loss 
and  did  not  involve  any  spurious  passenger  mutations  specific  to 
the  isolation  of  any  individual  nfl  mutant  line  (supplementary 
material  Fig.  S2).  Because  our  data  indicate  that  these  various  null 
alleles  are  equivalent,  we  refer  to  them  without  individual  allelic 
designations  henceforth  ( nfla” '”  and  nflb”1”). 

Mutants  carrying  at  least  one  wild- type  allele  of  either  nfla  or 
nflb  are  viable  and  fertile.  However,  when  crossing  parental 
genotypes  that  would  be  expected  to  yield  nfla”1”;  nflb progeny, 
none  were  observed  in  the  adult  population.  To  investigate  this 
further,  we  performed  quantitative  survival  studies.  At  7  dpf, 
nfla”1”;  nflb ”!”  larvae  began  to  die,  with  none  surviving  beyond 
10  dpf,  although  100%  of  wild-type  larvae  survived  to  10  dpf  (Fig. 
1G).  Furthermore,  100%  survival  at  10  dpf  was  also  observed  in 
nfla” 7~;  nflb*1*  («= 26),  nfla*1*;  nflb ”'”  («= 22),  nfla”1  ”;  nflb*1” 
(«= 28)  and  nfla*1”;  nflb”1”  («= 24)  larvae.  The  swim  bladders  of 
nfla”1”;  nflb”1”  larvae  were  frequently  observed  to  be  underinflated. 
However,  nfla”1”;  nflb”'”  larvae  maintained  the  ability  to  both 
consume  and  transit  live  paramecia,  suggesting  that  premature 
death  was  not  the  result  of  starvation  (supplementary  material  Fig. 
S3).  Additionally,  an  incompletely  penetrant  valvular  insufficiency 
phenotype  was  appreciated  in  nfla”1”;  nflb”1”  larvae,  as  well  as  in 
those  harboring  only  a  single  wild-type  nfl  allele  (nfla*1”;  nflb”1” 
and  nfla”1”;  nflb*1”)  (supplementary  material  Movies  1-3). 

OPC  and  Schwann  cell  hyperplasia  in  nfla”1”; nflb”1”  larvae 

We  previously  described  OPC  hyperplasia  after  nfla  and  nflb 
morpholino  knockdown  in  the  context  of  a  homozygous  pS3 
mutant  background  (Lee  et  al.,  2010).  To  examine  nfla  and  nflb 
function  in  OPCs  and  other  tissues  beyond  the  first  few  days  of 
life,  we  crossed  several  cell-type-specific  zebrafish  reporter  lines 
into  nfla/nflb  mutant  backgrounds.  At  2  dpf,  olig2  expression 
appeared  normal  in  Tg (olig2:GFP);  nfla”1”;  nflb”1”  embryos,  as 
assessed  by  both  whole-mount  in  situ  hybridization  analysis  of 
endogenous  olig2  mRNA  expression  and  GFP  expression 
(supplementary  material  Fig.  S4A-D).  We  also  evaluated  nfl  loss 
in  Tg (soxlO-.GFP)  embryos.  This  transgene  drives  GFP  expression 
in  specified  ventral  spinal  cord  OPCs,  but  not  the  neighboring 
motoneurons  that  arise  from  a  common  progenitor  cell,  as  well  as 
in  Schwann  cells  of  the  posterior  lateral  line  nerve  (PLLn).  At  2 
dpf,  similar  numbers  of  sox / 0: G F P- p o s i 1 1 ve  OPCs  were  detected 
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in  the  dorsal  and  ventral  spinal  cord  of  wild-type  and  nfla~ /_; 
nflb~  ~  embryos  (supplementary  material  Fig.  S4E,F).  Examination 
of  nfl- null  PNS  Schwann  cells  at  2  dpf  showed  no  effect  on  the 
number  of  .so,r/0:GFP-expressing  cells  associated  with  the  PLLn, 


which  innervate  skin  mechanosensory  neuromast  cells 
(supplementary  material  Fig.  S4G,H). 

However,  nfla nflb larvae  exhibited  increased  numbers  of 
OPCs  at  4  dpf  compared  with  controls,  as  evidenced  by  an  excess 
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Fig.  1.  ZFN  and  TILLING  strategies  generate  null  alleles  of  zebrafish  nfla  and  nflb.  (A,B)  Scheme  of  site  targeted  for  ZFN-mediated  cleavage  in  exon  26  of 
nfl  a  and  exon  17  of  nflb.  (C)  Alignments  of  nucleotide  sequences  from  wild-type  and  mutant  nfla  and  nflb  alleles.  ZFN  target  sites  appear  in  red.  Yellow  bases 
represent  insertions  that  arise  from  repair  of  ZFN-induced  mutagenic  lesions.  Black  bases  correspond  to  the  mutated  codon  in  the  nf1aLum  allele  generated  by 
a  TILLING  strategy,  with  the  mutated  base  appearing  in  blue.  (D)  ZFN-induced  mutagenic  lesions  in  nfla  and  nflb  produce  frameshift  mutations  that  lead  to 
truncated  protein  products  following  short  regions  of  altered  translation,  which  are  indicated  in  green.  The  nonsense  mutation  in  the  nf1aL1247X  allele  generated 
by  a  TILLING  strategy  also  appears  in  green.  (E)  The  truncated  protein  products  predicted  by  the  ZFN-  andTILLING-induced  mutant  nfla/nflb  alleles  all  harbor 
complete  or  partial  loss  of  the  neurofibromin  GAP-related  domain  (GRD).  (F)  Western  blot  analysis  for  neurofibromin  in  protein  lysates  from  3  dpf  wild-type, 
nflaA5/A5;nflb+/+,  nf1a+,+;  nfl fc>+10/+1°  and  nflaA5,A5;nflb+'0,+w  zebrafish  larvae  (lOOpg  each)  or  293T  cells  (25  pig)  demonstrates  absence  of  Nfl  protein  in  nfla~'~; 
nflb~'~  larvae.  Equal  loading  was  confirmed  by  stripping  the  membrane  and  reprobing  for  a-tubulin.  (G)  Kaplan-Meier  survival  analysis  demonstrates  that  100% 
of  nfl  a~'~;  nfl  b~'~  larvae  (n= 22)  die  by  10  dpf  as  compared  with  0%  of  wild-type  larvae  (n= 27). 

Disease  Models  &  Mechanisms  883 


Page  15 


Disease  Models  &  Mechanisms  •  DMM 


RESEARCH  ARTICLE 


Zebrafish  model  of  neurofibromatosis  type  1 


of  dorsally  migrated  o/zg2:GFP-positive  OPCs  (Fig.  2A,B)  along  with 
increased  numbers  of  both  dorsally  and  ventrally  positioned 
sox J0:G F P- p o si ti ve  OPCs  (Fig.  2C,D),  consistent  with  our  analyses 
of  nfl  morphants  at  3  dpf  (Lee  et  al.,  2010).  In  addition,  we  observed 
an  increase  in  soxl  0:GFP-positive  Schwann  cells  associated  with 
the  PLLn  (Fig.  2E,F).  This  increase  in  PLLn  Schwann  cell  number 
was  not  associated  with  altered  proliferation  of  these  cells 
(supplementary  material  Fig.  S5A-F).  To  assess  the  roles  of  nfla 
and  nflb  in  the  developing  radial  glial  cells  of  the  spinal  cord, 
nfla/nflb  mutants  were  crossed  into  the  Tg(gfap:GFP)  line,  which 
expresses  GFP  from  the  glial  fibrillary  acidic  protein  {gfap ) 
promoter.  At  4  dpf,  nfla~ nflb larvae  harboring  a  gfap'.GFP 
transgene  demonstrated  no  readily  discernible  differences  in 
,g/«/?:GFP- positive  spinal  cord  radial  glial  cells  as  compared  with 
wild-type  larvae  (supplementary  material  Fig.  S4I,J).  However,  gfap 
expression  in  the  extensive  processes  of  radial  glial  cells  precludes 
precise  quantification  and  might  obscure  subtle  differences. 

To  determine  whether  neuronal  numbers  increased  in  concert 
with  OPCs  in  nfla/nflb  mutant  larvae,  we  used  anti-HuC/D  and 
anti-SOXlO  antibodies  (see  Methods)  to  discriminate  between 
o%2:GFP-positive  neurons  and  OPCs,  respectively.  No  difference 
between  the  number  of  o %2 : G F P - p o s i t i ve  and  HuC/D-positive 
neurons  was  appreciable  in  4  dpf  spinal  cord  sections  from  wild- 
type  and  nfla~j~;  nflb larvae  (Fig.  2G,H,  green  and  magenta;  Fig. 
2K).  However,  the  numbers  of  o%2:GFP-positive  and  SoxlO- 
positive  OPCs  (Fig.  2G,H,  arrowheads;  Fig.  2L)  and  PLLn  Schwann 
cells  (Fig.  2I,J,M)  were  significantly  increased  at  4  dpf  in  nfla~j~; 
nflb larvae  relative  to  wild-type  controls.  OPC  cell  numbers 
continued  to  increase  in  nfla~!~;  nflb larvae  at  8  dpf,  as  reflected 
by  increased  numbers  of  dorsally  localizing  o  lig'2 :  G  F  P-  p  o  s  i  Li  ve 
OPCs  as  well  as  both  dorsally  and  ventrally  localizing  soxlO: GFP- 
positive  OPCs  (supplementary  material  Fig.  S4M-P).  Increased 
numbers  of  scwIO:GFP-expressing  PLLn  Schwann  cells  were  also 
evident  in  nfla~j~;  nflb larvae  at  8  dpf  (supplementary  material 
Fig.  S4Q,R).  Thus,  loss  of  nfla  and  nflb  does  not  affect  the 
specification  of  OPCs  at  2  dpf,  but  instead  promotes  the  progressive 
expansion  of  OPCs  without  a  concomitant  increase  in  neuronal 
cell  numbers.  Furthermore,  nfla/nflb  loss  triggers  Schwann  cell 
hyperplasia  beginning  at  4  dpf. 

Immunohistochemical  analysis  using  the  Zrfl  antibody,  which 
labels  Gfap  in  zebrafish,  showed  coexpression  with  GFP  expressed 
from  the  gfap:GFP  transgene  and  revealed  a  similar  pattern  of 
expression  in  wild-type  and  nfla~'~;  nflb~'~  larvae  at  4  dpf 
(supplementary  material  Fig.  S6A,B).  Zebrafish  radial  glial  cells  also 
express  brain  lipid-binding  protein  (Bibp),  and 
immunohistochemical  analysis  of  wild-type  and  nfla~j~;  nfl b^1 
spinal  cords  with  an  anti-BLBP  antibody  at  4  dpf  revealed  an 
obvious  decrease  in  Blbp  expression  in  the  ^n/JiGFP- positive  radial 
glia  of  nfla nflb larvae  (supplementary  material  Fig.  S6C,D). 
These  results  suggest  that  although  g/B/>:GFP-positive  radial  glial 
cells  in  nflb larvae  appear  normal  in  number,  they  fail 

to  express  appropriate  levels  of  Blbp  indicating  a  defect  in 
gliogenesis.  An  additional  abnormality  of  gliogenesis  was  observed 
at  8  dpf  as  a  disruption  in  the  regular  segmental  pattern  of  glial 
process  outgrowth  in  Tg {gfap:GFP);  nfla~'~;  nflb~'~  larvae 
(supplementary  material  Fig.  S4K,L). 
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Fig.  2.  Loss  of  nfl  a  and  nflb  causes  hyperplasia  of  OPCs  and  Schwann  cells. 

(A,B)  Confocal  images  of  spinal  cords  in  nf1a~'~;  nfl  b~';Tg(olig2:GFP)  larvae  (B) 
demonstrate  increased  numbers  of  dorsally  migrating  (white  brackets) 
o//g2:GFP-positive  OPCs  as  compared  with  wild-type;  Tg(olig2:GFP)  larvae  (A)  at  4 
dpf.  (C,D)  Confocal  images  of  spinal  cord  in  nfl a~'~;  nfl b~'~;Tg(sox10:GFP )  larvae 
(D)  demonstrate  increased  numbers  of  both  dorsally  (white  brackets)  and 
ventrally  (red  brackets)  positioned  soxl 0:GFP-positive  OPCs  as  compared  with 
wild-type;  Tg  (sox  7  O.-GFP)  larvae  (C)  at  4  dpf.  (E,F)  nfl  a~'~;  nfl  tf/_;Tg(sox7  0:GFP) 
larvae  (F)  show  an  increased  number  of  sox70:GFP-positive  Schwann  cells 
associated  with  the  peripheral  lateral  line  nerve  (PLLn;  white  brackets)  as 
compared  with  wild-type;  Jg(sox10:GFP)  larvae  (E)  at  4  dpf.  (G,H)  Neuronal 
numbers  (o//g2:GFP-  [green],  FluC/D-  [magenta]  positive)  do  not  increase  in 
concert  with  OPCs  (olig2: GF-  [green],  Soxl  0-  [red]  positive;  arrowheads)  in 
transverse  sections  through  the  spinal  cord  of  nfl  or1';  nfl  b~/~;Tg{olig2:GFP) 
larvae  (H)  as  compared  to  wild-type;  Tg (olig2:GFP)  larvae  (G)  at  4  dpf. 

(I,J)  Increased  numbers  of  PLLn  Schwann  cells  (sox70:GFP-  [green],  SoxlO-  [red] 
positive)  are  appreciated  in  transverse  sections  of  nfl  a'~;  nfl  tf/_;Tg(sox70:GFP) 
larvae  (J)  as  compared  with  wild-type;  Jg(sox10:GFP)  (I)  larvae  at  4  dpf. 

(K,L)  Quantification  of  neurons  (o//g2:GFP-,  FluC/D-positive  cells)  (K)  and  OPCs 
(o//'g2:GFP-,  Soxl  0-positive  cells)  (L)  from  transverse  sections  through  the  spinal 
cord  of  wild-type;  ~fg(olig2:GFP)  and  nf1a~,~;nf1b~,~;Jg(olig2:GFP)  larvae  at4dpf. 
Values  indicate  mean  +  s.e.m.  per  section  (n=30  from  five  each  of  wild-type  and 
nfl cf1'; nf1b~'~  larvae).  (M)  Quantification  of  sox70:GFP-positive  Schwann  cells  in 
the  PLLn  of  wild-type;  Tg  (sox  7  O.GFP)  and  nflcr'-tnflF^^gisoxIOtGFP)  larvae  at 
4  dpf.  Values  indicate  mean  +  s.e.m.  per  hemisegment  (n= 5  each  for  wild-type 
and  nf1a~'~;  nf1b~'~  larvae).  ***P<0.001 .  Scale  bars:  20  pm. 
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OPC  hyperplasia  in  nfla~h;  nflb~'~  larvae  results  from  increased 
proliferation 

In  njla1;  nflb  larvae,  OPC  numbers  were  increased  relative  to 
control  animals  at  4  dpf,  but  not  at  2  dpf,  suggesting  that  nfl- null 
OPCs  might  proliferate  faster  during  this  time  period.  We  assessed 
OPC  proliferation  by  labeling  larvae  with  BrdU  for  12  hours  starting 
at  3.5,  4.5  and  5.5  dpf  followed  by  immunohistochemical  analysis 
of  transverse  sections  through  the  spinal  cord.  By  4  dpf  and 
continuing  through  6  dpf,  the  number  of  o%2:GFP-positive  and 
SoxlO-positive  OPCs  in  nfl- null  larvae  was  significantly  increased 
in  comparison  with  wild-type  controls  (Fig.  3A-H,  arrows;  Fig.  31). 
Mutant  larvae  further  exhibited  significantly  increased  numbers 
of  olig2: GFP-,  SoxlO-  and  BrdU-positive  OPCs  at  4  and  5  dpf  as 
compared  with  controls,  indicating  that  neurofibromin  normally 
suppresses  the  proliferation  of  OPCs  during  this  period  of 
development  (Fig.  3A-H,  arrowheads;  Fig.  3J).  There  was  little 
detectable  BrdU  incorporation  at  6  dpf  in  either  mutant  or  control 
populations. 

Myelination  is  aberrant  in  nfl  a1;  nflb~'~  larvae 

We  evaluated  the  ability  of  nfla~!~;  njlb  OPCs  to  differentiate 
appropriately  by  examining  the  gene  expression  levels  of 
proteolipid  protein  la  (plpla )  and  myelin  basic  protein  ( mbp ), 
markers  of  differentiated  oligodendrocytes,  in  wild-type  and 
mutant  larvae  at  5  dpf.  In  nfl- null  larvae,  fewer  plpla- positive 
cells  were  detected  in  the  midbrain  and  hindbrain  regions  (Fig. 
4A,B)  as  well  as  along  the  dorsal  and  ventral  spinal  cords  (Fig. 
4C,D)  as  compared  with  controls.  Central  nervous  system  (CNS) 


Fig.  3.  Increased  proliferation  drives  OPC  hyperplasia  in  nf1a~'~; nflb~'~ 
larvae.  (A-H)  Transverse  spinal  cord  sections  of  wild-type;  Jg(olig2:GFP)  (A-D) 
and  nfl a'~; nfl b~/~;Tg(olig2:GFP)  larvae  (E-H)  labeled  with  anti-BrdU  antibody 
(B,F,  red)  or  anti-SOXIO  antibody  (C,G,  magenta)  at  4  dpf.  Arrows  indicate 
BrdU-negative,Tg(o//g2:GFP)-,  SoxlO-positive  OPCs.  Arrowheads  indicate 
BrdU-positive,  Tg (olig2:GFP)-,  Soxl  0-positive  OPCs.  (I,J)  Quantification  of  total 
(I)  and  BrdU-positive  OPCs  (J)  from  transverse  spinal  cord  sections  of  wild- 
type;  Tg  (olig2:GFP)  and  nf1a~,~;nf1b~/  ;Jg(olig2:GFP)  larvae  at  4,  5  and  6  dpf. 
Values  indicate  mean  +  s.e.m.  per  section  (n= 30  from  five  each  of  wild-type 
and  nffo_/“;nfft>_/“  larvae).  **P<0.01;***P<0.001.  Scale  bars:  20  pm. 

Disease  Models  &  Mechanisms 


expression  of  mbp,  on  the  other  hand,  was  indistinguishable 
between  wild-type  and  nfla~j~;  nflb larvae  at  5  dpf 
(supplementary  material  Fig.  S7).  However,  mbp  expression  was 
elevated  in  Schwann  cells  of  the  PLLn  in  nfla~'~;  nflb  larvae 
as  compared  with  controls  (Fig.  4E,F,  arrowheads).  These  data 
are  consistent  with  perturbed  oligodendrocyte  differentiation  in 
the  CNS  as  well  as  in  PNS  Schwann  cells  associated  with  the  PLLn 
of  nfla~ nflb larvae. 

We  went  on  to  examine  the  ultrastructure  of  myelinated  CNS 
axons  in  control  and  nf  1-null  larvae  by  transmission  electron 
microscopy  (TEM).  At  8  dpf,  oligodendritic  myelin  sheaths  were 
tightly  wrapped  around  CNS  axons  in  the  ventral  spinal  cord  of 
wild-type  larvae  (Fig.  4G,I).  By  contrast,  nflb axons  were 

loosely  encircled  by  multiple  lamellae  rather  than  by  compact 
myelin  sheaths,  indicating  that  neurofibromin  is  required  for  the 
normal  formation  of  the  concentric  layers  of  oligodendrocyte 
membranes  that  enwrap  neuronal  axons  of  the  CNS  to  promote 
neural  conduction  (Fig.  4H,J). 

wild  type  nflar'  -,  nflb v' 


5  dpf 


Fig.  4.  nfl o+; nfl b~'~  larvae  exhibit  myelination  defects.  (A-D)  plpla 
expression  is  decreased  in  glial  cells  of  the  midbrain  and  hindbrain  regions 
(A,B,  boxes)  as  well  as  the  dorsal  (black  bracket)  and  ventral  (white  brackets) 
spinal  cord  of  nfl  a~'~;  nfl  b~'~  larvae  (B,D)  as  compared  to  wild-type  larvae  (A,C) 
by  whole-mount  in  situ  hybridization  at  5  dpf.  (E,F)  mbp  expression  is  elevated 
in  Schwann  cells  of  the  PLLn  (arrowheads)  of  nfl  a~'~;  nfl  b~'~  larvae  (F)  as 
compared  to  wild-type  larvae  (E)  by  whole-mount  in  situ  hybridization  at  5 
dpf.  (G-J)  Transverse  sections  through  the  ventral  spinal  cord  nf  1  a~f~;  nf  1  b~'~ 
larvae  (H,  boxed  region  magnified  in  J)  demonstrate  defects  in  formation  of 
compact  myelin  sheaths  as  compared  with  wild-type  larvae  (G,  boxed  region 
magnified  in  I)  at  8  dpf.  Scale  bars:  1 00  pm  (A-F)  and  0.5  pm  (G-J). 
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Loss  of  nfl a  and  nflb  causes  upregulation  of  Ras  signaling  in  the 
spinal  cord 

Given  the  well-described  role  of  neurofibromin  as  a  negative 
regulator  of  Ras,  we  hypothesized  that  nfl  loss  in  our  mutants  would 
lead  to  activation  of  downstream  effector  pathways.  Western  blot 
analysis  of  whole  larvae  extracts  revealed  an  upregulation  of 
phosphorylated  ERK1  and  ERK2  (pERKl/2)  in  nfla- nflb -l- 
larvae  at  3  dpf,  whereas  levels  in  nflaT1-;  nflb+,+  and  nfla+,+;  nflb-1- 
larvae  remained  unchanged  (Fig.  5A).  These  data  are  consistent 
with  the  absence  of  functioning  neurofibromin  protein  in  nfla-,~; 
nflb-1-  larvae  and  support  functional  redundancy  between  nfl  a 
and  nflb.  We  next  assessed  the  activation  of  Ras  effector  pathways 
in  the  spinal  cords  of  wild-type,  nfla+,~;  nflb-1-,  nfla-1-;  nflb*1- 
and  nfla-1-;  nflb-1-  animals  by  immunohistochemical  analysis  of 
transverse  larval  sections.  Antibodies  directed  against  HuC/D, 
pERI<l/2,  and  phosphorylated  S6  (pS6)  were  used  to  label  neurons 
and  assess  activation  of  ERK  and  mTOR  signaling  pathways, 
respectively  (Fig.  5B-Q).  Although  pERKl/2  staining  was  only 
minimally  observed  in  a  few  neurons  and  portions  of  spinal  cord 
white  matter  at  4  dpf  in  wild-type  larvae  (Fig.  5B,N),  a  striking 
upregulation  of  pERKl/2  was  detected  in  nfla-1-;  nflb-1-  larvae 
(Fig.  5E,Q).  Increased  ERK  signaling  was  also  noted  at  3  dpf  in  spinal 
cord  neurons  and  white  matter  of  nfla-1-;  nflb-1-  larvae,  but  was 
absent  at  2  dpf  (supplementary  material  Fig.  S8).  Although  pS6 
signaling  was  evident  in  multiple  spinal  cord  neurons,  we  observed 
no  differences  in  these  cells  between  wild-type  and  nfl  mutant 
animals  at  2,  3  or  4  dpf  (supplementary  material  Fig.  S8Ae-Ah,Be- 
Bh;  Fig.  5F-I).  These  data  suggest  that  activation  of  mTOR  signaling 
(as  assessed  by  S6  phosphorylation)  is  not  altered,  at  least  in  the 
spinal  cord  following  nfl  loss. 


nfl  and  p53  cooperate  to  accelerate  zebrafish  tumorigenesis  in 
vivo 

Mammalian  NF1  has  been  shown  to  be  a  potent  tumor  suppressor; 
however,  we  did  not  identify  any  tumors  over  18  months  of 
observation  in  adult  zebrafish  homozygous  for  either  mutant  nfl 
allele  alone  ( nfla nflb*1*  or  nfla*1*;  nflb-1-)  or  in  combination 
with  heterozygous  loss  of  the  remaining  allele  ( nfla *1-;  nflb-1-  or 
nfla-1-;  nflb*1-).  Loss  of  p53  has  been  shown  to  cooperate  with 
NF1  (Cichowski  et  al.,  1999;  Vogel  et  al.,  1999)  as  well  as  other 
mutations  that  activate  Ras  signaling  (Eliyahu  et  al.,  1984;  Parada 
et  al.,  1984;  Kemp  et  al.,  1994;  Tanaka  et  al.,  1994;  Hundley  et  al., 
1997)  in  mammalian  tumorigenesis,  so  we  next  bred  pS3  mutant 
zebrafish  into  an  nfl- mutant  background  to  generate  nfla*1-; 
nflb-1-;  p53e7le7  fish.  These  animals  were  incrossed  to  derive 
nfla*'*;  nflb-1-;  p53e7le7  and  nfla*1-;  nflb-1-;  pS3e7le7  fish,  which 
were  subsequently  monitored  carefully  for  tumorigenesis.  At  31 
weeks  post  fertilization  (wpf),  nfla*1-;  nflb-1-;  p53e7le7  fish  began 
to  develop  tumors  with  high  penetrance  although  only  one  nfla*1*; 
nflb-1-;  p53e7le7  fish  developed  a  tumor  at  44  wpf  (Fig.  6A).  At  45 
wpf,  tumor  penetrance  was  higher  in  nfla*1-;  nflb-1-;  p53e7le 7  fish 
(24/39;  62%)  than  in  nfla*'*;  nflb-1-;  p53e7le7  fish  (1/14;  7%).  We 
have  previously  reported  that p53e7le7  fish  with  wild-type  nfl  alleles 
develop  MPNSTs.  These  tumors  did  not  begin  to  develop  until  41 
weeks  of  age,  which  was  similar  to  the  results  with  nfla*1*;  nflb-1-; 
p53e7,e7  animals  (Fig.  6A).  Furthermore,  the  penetrance  of  tumors 
in  p53-null  animals  was  only  28%  by  66  wpf  (Berghmans  et  al., 
2005).  Thus,  the  combined  loss  of  p53  and  3  of  4  nfl  alleles  in 
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Fig.  S.nflalnflb  mutants  exhibit  upregulation  of  pERK1/2.  (A)  Western 
blot  analysis  for  pERK1/2  in  protein  lysates  prepared  from  wild-type,  nfla 
nfl  b*/+,  nfl  a*'*;  nfl  b-'-  and  nfla-1-;  nflb-/-  larvae  (1 00  pig)  reveals  increased 
pERK1/2  levels  in  nfla-/-;nflb~/ ~  larvae  as  compared  with  wild-type,  nfla-'-; 
nflb+/+  and  nfla+,+;nflb-/-  larvae  at  3  dpf.  Equal  loading  was  confirmed  by 
stripping  the  membrane  and  reprobing  for  total  ERK1/2.  (B-Q)  Transverse 
spinal  cord  sections  of  wild-type  (B,F,J,N),  nfl  a*/-;  nflb-'-  (C,G,K,0),  nfla-1-; 
nflb+/-( D,H,L,P)  and  nfla-1-;  nflb-'-  (E,I,M,Q)  larvae  labeled  with  anti-pERK1/2 
antibody  (B-E,  green),  anti-pS6  antibody  (F-l,  red),  or  anti-HuC/D  antibody  (J-M, 
magenta)  demonstrate  marked  upregulation  of  pERK1/2  in  nfla-'-; nfl b-1- 
larvae  and  intermediate  levels  of  pERK1/2  in  larvae  harboring  a  single 
functioning  nfl  allele  as  compared  with  wild-type  larvae  at  4  dpf  (n=5  each  for 
wild-type  and  mutant  larvae).  Scale  bars:  40  pm. 


zebrafish  markedly  accelerates  the  onset  and  increases  the 
penetrance  of  tumors  as  compared  with  the  loss  of  p53  alone  or 
the  concomitant  loss  of  p53  and  both  alleles  of  nflb,  but  with  intact 
nfla. 

Tumors  in  nfla*1-;  nflb-1-;  p53e7le7  fish  were  observed  in  the 
brain  («= 2),  eye  («= 8),  gill  («=  1),  abdomen  («= 8)  and  trunk  («= 5). 
Brain  tumors  developed  very  early  (3 1  and  33  wpf)  (Fig.  6A,  arrows) 
and  demonstrated  features  of  diffuse  high-grade  gliomas,  whereas 
all  other  tumor  types  were  most  consistent  with  MPNSTs  (Fig.  6B- 
G).  Histopathologically,  the  brain  tumors  were  highly  cellular  and 
composed  of  ovoid  to  rounded  cells  with  marked  nuclear 
pleomorphism  and  diffuse  single  cell  infiltration  of  parenchyma, 
including  pre-existing  neurons  (Fig.  7A-F).  Occasional  mitoses  were 
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nfla w-;  nflb';  p53P7/e7 

glioma  (31  weeks)  MPNST  (32  weeks) 


Fig.  6.  nflalnflb  mutants  demonstrate  increased  susceptibility  to 
tumorigenesis  in  a  pS3  mutant  background.  (A)  Kaplan-Meier  tumor-free 
survival  analysis  for  nf'la+l~;nf'lb~l~;p53e:7l<! 7  (n= 39)  animals  demonstrates 
significantly  decreased  survival  as  compared  with  nfl a+,+;  nfl b~'~; p53e7,e7 
(n=14)  animals  (P<0.001).  Arrows  indicate  ages  at  which  animals  were 
identified  with  brain  tumors  demonstrating  features  of  diffuse  high-grade 
gliomas.  (B,C)  31-week-old  nfla+,~;nf1b~,~;p53e7/!!7 animal  with  a  high-grade 
glioma  (B,  arrowhead)  and  32-week-old  nfl  a+,~;  nfl  b^'~;  p53e7,e7  anima\  with  a 
malignant  peripheral  nerve  sheath  tumor  (MPNST)  (C,  arrowhead).  (D-G)  H&E 
staining  of  sagittal  sections  through  the  high-grade  optic  glioma  (D,  boxed 
area  magnified  in  F)  or  MPNST  (E,  boxed  area  magnified  in  G).  Scale  bars: 

200  pm  (D,E)  and  50  pm  (F,G). 


identified,  but  no  necrosis  or  vascular  proliferation  was  detected 
(Fig.  7G-I).  Assessment  of  tumor  lineage  by  immunohistochemical 
analysis  showed  that  approximately  80%  of  tumor  cells  stained 
positive  for  the  oligodendroglial  marker  SoxlO,  with  little  staining 
in  matched  wild-type  tissue  (Fig.  7J-L).  The  presence  of  a  SoxlO- 
negative  tumor  cell  subpopulation  is  consistent  with  the  level  of 
heterogeneity  for  oligodendroglial  transcription  factors,  such  as 
SoxlO  and  01ig2,  and  is  characteristic  of  astrocytic  or  mixed  gliomas 
as  compared  with  pure  oligodendroglial  class  tumors  (Fig.  7K, 
arrowhead)  (Ligon  et  al.,  2004;  Bannykh  et  ah,  2006).  Staining  for 
the  astrocytic  marker  Gfap  (Fig.  7M-0)  highlighted  a  subpopulation 
of  cells  within  the  tumor  with  coarse,  irregular  cytoplasmic 


Fig.  7.  Tumors  from  nfl a+l~;  nfl b~'~; pS3s7le7  animals  express  glial  markers. 

(A-F)  H&E  staining  of  sagittal  sections  through  p53e7/e7  brain  tissue  (A,D)  or 
brain  tumors  in  nf'la+l~;nf1b~l~;p53^7le7 a nimals  at  31  (B,E)  and  33  (C,F)  wpf 
(cb,  cerebellum;  v,  ventricle).  The  arrowhead  in  E  shows  infiltration  of  a  single 
cell  through  the  parenchyma  and  around  normal  neurons. 

(G-U)  Immunohistochemical  analyses  for  phosphorylated  histone  H3  (G-l, 
pH3),  SRY-box  10  (J-L,  SoxlO),  glial  fibrillary  acidic  protein  (M-O,  Gfap),  HuC/D 
(P-R,  Hu)  and  synaptophysin  (S-U,  Sypa)  on  sagittal  sections  through  p53e7,e7 
brain  tissue  (G,J,M,P,S)  or  brain  tumors  in  nf1a+,~; nf1b~'~; p53e7,e7  animals  at  31 
(H,K,N,Q,T)  and  33  (l,L,0,R,U)  wpf  demonstrate  occasional  mitoses  in  tumor 
tissue  (H,l)  with  most  tumor  cells  staining  positive  for  the  oligodendroglial 
marker  SoxlO  (K,L),  a  Soxl 0-negative  subpopulation  (K,  arrowhead),  an 
astrocytic  component  (N,  arrowheads  and  O)  and  the  absence  of  mature 
neuronal  markers  HuC/D  (Q,R)  or  synaptophysin  (T,U).  Scale  bars:  450  pm  (A-C) 
and  20  pm  (D-U). 


processes  also  consistent  with  the  presence  of  an  astrocytic  lineage 
component  (Fig.  7N,  arrowheads).  Tumor  cells  did  not  express  the 
mature  neuronal  markers  HuC/D  or  synaptophysin  (Sypa), 
consistent  with  their  glial  origin  (Fig.  7P-U). 
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Immunohistochemical  analysis  for  pERI<l/2  and  pS6  to  assess 
activation  of  ERK  and  mTOR  signaling  pathways,  respectively 
(supplementary  material  Fig.  S9A-F)  revealed  increased  pERI<l/2 
staining  in  the  more  malignant  of  these  two  brain  tumors,  with 
normal  amounts  of  pS6  (supplementary  material  Fig.  S9C,D).  The 
more  hyperplastic  brain  tumor  showed  increased  pERKl/2  and  pS6 
staining  (supplementary  material  Fig.  S9E,F).  These  data 
demonstrate  that  some,  but  not  all,  brain  tumors  in  nfla +l-;  nflb-1-; 
p53e7,e7  animals  demonstrate  hyperactivation  of  ERK  and  mTOR 
pathways,  consistent  with  mouse  and  human  AIFl-derived  MPNSTs 
and  gliomas  (Dasgupta  et  al.,  2005;  Zhu  et  al.,  2005b).  Collectively, 
these  findings  suggest  that  the  tumors  were  high-grade  gliomas 
most  closely  resembling  human  anaplastic  astrocytoma  or 
anaplastic  oligoastrocytoma  WHO  grade  III.  Furthermore, 
MPNSTs  (Fig.  6C,E,G)  exhibited  spindle-shaped  tumor  cells  and 
extensive  necrosis  consistent  with  this  tumor  type  (Ducatman  et 
al.,  1986;  Wanebo  et  al.,  1993;  Hirose  et  al.,  1998).  Taken  together, 
we  conclude  that  nfla  and  nflb  mutations  cooperate  with p53  loss 
to  generate  high-grade  gliomas  and  MPNSTs. 

nfla-1-;  nflb-1-  larvae  show  motor  and  learning  deficits 

Deficits  in  motor  coordination  and  cognition,  including  learning 
and  memory,  are  characteristic  of  NF1  patients  and  animal  models. 
To  examine  motor  behavior  and  cognition  in  nfla -l-;  nflb larvae, 
we  performed  kinematic  analysis  of  the  short-latency  C-start 
(SLC),  a  highly  stereotyped  yet  modifiable  acoustic  startle  reflex  in 
the  zebrafish  (Burgess  and  Granato,  2007a;  Wolman  et  al.,  2011). 
Unlike  their  siblings,  nfla-1-;  nflb-1'  larvae  showed  a  deficit  in 
short-term  SLC  habituation  when  presented  with  repetitive 
acoustic  stimulation  (supplementary  material  Fig.  S10A). 
Furthermore,  nfla nflb-1-  larvae  performed  kinematically 
weaker  SLC  responses,  as  indicated  by  decreased  head  turning 
angle,  maximum  angular  velocity  and  distance  traveled  following 
delivery  of  an  acoustic  stimulus  (supplementary  material  Fig.  S10B- 
D).  Taken  together,  these  data  support  the  hypothesis  that  nfla  and 
nflb  function  redundantly  during  zebrafish  development  and  that 
only  one  of  the  four  wild-type  nfl  alleles  is  required  for 
phenotypically  normal  embryonic  development,  motor  behavior 
and  cognition. 

nfl  mutants  exhibit  melanophore  defects 

Notably,  nfl  null  larvae  displayed  aberrant  lateral  stripe 
pigmentation  as  compared  with  wild-type  controls  at  6  dpf  (Fig. 
8A-D).  This  phenotype  was  first  appreciable  at  4  dpf  and  was 
manifested  as  a  disruption  in  the  uniform  pattern  of  melanophores 
arranged  along  the  lateral  stripes  (Fig.  8B,D,  brackets).  To  further 
investigate  this  phenotype,  we  quantified  melanophore  numbers 
along  the  lateral  stripes  of  3  and  6  dpf  wild- type  and  nfla -/~;  nflb-1- 
larvae.  No  significant  differences  in  the  number  of  melanophores 
comprising  the  3  dpf  lateral  stripes  were  appreciable  in  any 
combination  of  mutant  nfl  alleles  as  compared  with  sibling  wild- 
type  controls  (supplementary  material  Fig.  Sll).  However,  at  6  dpf, 
nfla nflb-1-  larvae  exhibited  a  significant  reduction  in  the 
number  of  lateral  stripe  melanophores  (Fig.  8E).  Less  severe,  but 
still  significant,  decreases  were  also  noted  in  larvae  carrying  two 
or  three  mutant  nfl  alleles  (Fig.  8E).  No  difference  in  the  number 
of  apoptotic  cells  between  wild-type  and  mutant  larvae  was 
discernible  at  3  dpf  (68.1±11  cells/larva,  «=  12  wild-type  versus 
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Fig.  8.  nfl  mutants  display  decreased  lateral  stripe  melanophores. 

(A-D)  Lateral  stripe  melanophores  of  wild-type  (A,B)  larvae  demonstrate  a 
normal  uniform  pattern  that  is  absent  in  nfla-'-;  nfl  b-,-\ar\/ae  (C,D)  at  6  dpf. 
The  brackets  in  D  highlight  regions  where  melanophores  are  absent. 

(E)  Quantification  of  lateral  stripe  melanophores  from  wild-type  (n= 26), 
nfla+,~;  nflb-/-  (n= 30),  nfla-/-;nf1b+,+  (n=14),  nfla-,-;nflb+,~  (n= 26)  and  nfla 
nflb (n=12)  larvae  at  6  dpf.  Each  point  represents  the  number  of  lateral 
stripe  melanophores  in  an  individual  embryo  and  red  lines  indicate  mean 
values.  *P< 0.05,  **P< 0.01 .  Scale  bars:  0.5  mm  (A,C)  and  1 50  pm  (B,D). 


63.8±8.5,  «=12  mutants;  P= 0.27)  (supplementary  material  Fig. 
S12A-F).  Serial  observation  of  melanophore  development  in 
individual  nfla-1-;  nflb-1-  and  wild- type  larvae  revealed  a  defect 
of  migration  or  differentiation  of  regeneration  and  metamorphic 
lineage  melanophores.  At  3  dpf,  a  time  point  at  which  embryonic 
melanophore  development  is  complete  (Hultman  et  al.,  2009),  no 
difference  in  melanophore  number  along  the  lateral  stripe  was 
observed  (supplementary  material  Fig.  S13A,B,G),  indicating  that 
the  melanophores  develop  normally  in  nfl  mutants.  To  assess  the 
regeneration  and  metamorphic  lineage  of  melanophores,  we 
suppressed  melanin  synthesis  after  3  dpf  by  treatment  with  N- 
phenylthiourea  (PTU)  (Hultman  and  Johnson,  2010),  which  allowed 
us  to  identify  newly  formed  melanophores  by  their  melanin¬ 
negative  and  pale  appearance  (supplementary  material  Fig.  S13C,D, 
red  arrows).  Removal  of  PTU  at  5  dpf  restored  melanin  synthesis 
and  regeneration  and  metamorphic  melanophores  appeared 
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melanin-positive  by  6  dpf  (supplementary  material  Fig.  S13E,F, 
black  arrows).  The  abnormal  appearance  of  the  lateral  stripes  in 
mutant  larvae  can  be  attributed  to  defects  in  patterning  of 
regeneration  and  metamorphic  melanophores  (supplementary 
material  Fig.  S13F,G,  brackets),  suggesting  abnormal  migration  or 
differentiation  of  this  lineage.  Head  melanophore  numbers  at  6  dpf 
were  not  significantly  changed  in  nfl  mutant  larvae  and  sibling 
controls  (supplementary  material  Fig.  S14).  Collectively,  these  data 
demonstrate  a  specific  defect  in  lateral  stripe  melanophore  numbers 
following  nfl  allele  loss,  most  prominent  in  the  setting  of  biallelic 
nflalnflb  loss.  Each  nfla~!~;  nflb larvae  showed  a  unique 
pattern  of  lateral  stripe  melanophore  loss,  suggesting  a  stochastic 
defect  in  melanophore  differentiation  from  neural  crest, 
proliferation,  migration  or  survival,  rather  than  a  defect  specific  to 
any  particular  somite  of  the  developing  embryo  and  larvae. 

DISCUSSION 

NF1  is  a  common  tumor- predisposing,  autosomal  dominant  genetic 
disorder  characterized  by  cafe-au-lait  macules  and  cutaneous 
neurofibromas.  In  addition  to  these  pigmentation  defects  and 
tumors  of  the  peripheral  nervous  system  (PNS),  NF1  patients 
demonstrate  highly  diversified  clinical  features  with  multiple  tissue 
types  affected.  Prominent  lesions  include  optic  pathway  gliomas, 
Lisch  nodules,  skeletal  dysplasia,  cardiovascular  abnormalities, 
learning  defects  and  various  cancers  such  as  leukemia  and  intestinal 
tumors  (Side  et  al.,  1997;  Side  et  al.,  1998;  Bahuau  et  al.,  2000; 
Andersson  et  ah,  2005;  Williams  et  ah,  2009).  It  remains  unclear 
how  NF1  mutations  contribute  to  the  diverse  symptoms  and  tissue 
types  affected  in  patients.  Because  neurofibromin  is  a  very  large 
protein  that  is  highly  conserved  evolutionarily,  it  is  likely  to  have 
activities  related  to  functional  domains  other  than  those  affecting 
Ras  signaling.  For  example,  in  addition  to  Ras,  neurofibromin  can 
bind  to  microtubules,  syndecan,  phospholipids  and  amyloid 
precursor  protein  (Xu  and  Gutmann,  1997;  Hsueh  et  ah,  2001; 
D’Angelo  et  ah,  2006;  De  Schepper  et  ah,  2006).  Intriguingly,  recent 
studies  indicate  that  neurofibromin  might  function  as  a  positive 
regulator  of  adenylyl  cyclase  (Guo  et  ah,  2000;  Tong  et  ah,  2002; 
Dasgupta  et  ah,  2003).  This  function  of  neurofibromin  modulates 
neuronal  differentiation  (Hegedus  et  ah,  2007),  suggesting  the 
possibility  that  cognitive  defects  in  NF1  patients  might  be  related 
to  defects  in  cAMP  signaling  rather  than  activated  Ras. 
Furthermore,  the  tremendous  clinical  variability  in  the  phenotypic 
spectrum  seen  among  families  with  the  same  molecular  NF1  lesion 
posits  the  role  of  unlinked  modifier  loci  in  regulating  the 
expressivity  of  disease  characteristics  (Easton  et  ah,  1993;  Sabbagh 
et  ah,  2009).  Nevertheless,  the  identification  of  specific  modifier 
genes  and  the  relative  contributions  of  Ras  signaling  versus  other 
neurofibromin-regulated  pathways  for  specific  phenotypes  have  yet 
to  be  fully  elucidated.  The  zebrafish  model  that  we  have  developed 
offers  an  attractive  tool  for  furthering  this  analysis  because  it  is 
amenable  to  small  molecule  screens,  genetic  modifier  screens  and 
genetic  rescue  experiments. 

A  hallmark  feature  of  human  NF1  is  the  presence  of  pigmentation 
defects  known  as  cafe-au-lait  spots.  Little  is  known  about  the 
underlying  mechanisms  responsible  for  this  abnormality,  and 
pigmentation  abnormalities  in  other  animal  models  of  NF1  have 
not  been  described.  In  this  regard,  it  is  of  interest  that  zebrafish 
lacking  neurofibromin  exhibit  abnormal  patterning  of  the 
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melanophores  that  compose  the  lateral  stripes,  a  phenotype  similar 
to  that  following  pharmacologic  inhibition  of  the  upstream  Ras 
effector  ErbB  (Hultman  et  al.,  2009;  Hultman  and  Johnson,  2010). 
This  easily  observable  and  quantifiable  phenotype  offers  the 
opportunity  to  probe  underlying  molecular  pathways  modulated 
by  nfl  in  melanophores. 

Several  studies  employing  murine  models  have  previously  shown 
a  role  for  neurofibromin  in  regulating  cell  numbers  of  CNS  OPCs 
and  astrocytes  as  well  as  PNS  Schwann  cells  and  sympathetic 
neurons  (Brannan  et  al.,  1994;  Gutmann  et  al.,  1999;  Bennett  et  al., 
2003;  Zhu  et  al.,  2005b;  Hegedus  et  al.,  2007;  Zheng  et  al.,  2008). 
These  are  consistent  with  our  findings  in  nfla/nflb  mutants  and 
indicative  of  a  strong  evolutionary  conservation  of  Nfl  function  in 
neural  development.  Our  observation  of  impaired  compact  myelin 
formation  and  reduced  CNS  plpla  expression  with  unperturbed 
CNS  and  increased  PNS  mbp  expression  suggests  that  the 
differentiation  programs  of  oligodendrocytes  and  Schwann  cells 
respond  differently  to  Nfl  deficiency.  Alternatively,  the  observed 
ultrastructural  defects  in  myelinated  CNS  axons  might  arise 
secondary  to  neuronal  defects.  The  accessibility  of  the  zebrafish 
embryo  to  mosaic  analysis  offers  the  ability  to  differentiate  between 
these  possibilities.  Unlike  Nfl -deficient  mice,  however,  nflalnflb 
mutant  radial  glia  failed  to  demonstrate  an  appreciable  increase  in 
Gfap  expression.  Instead,  we  observed  a  decrease  in  Blbp  expression 
and  irregularities  in  patterning  of  nfla~'~;  nflb ~/_  radial  glial  cells. 
This  discrepancy  might  reflect  species-specific  differences  in  neural 
tissue  as  opposed  to  Nfl  function,  because  it  remains  unclear 
whether  zebrafish  radial  glia-like  ependymal  cells  are  functionally 
equivalent  to  mammalian  astrocytes. 

We  have  previously  characterized  cardiovascular  defects 
resulting  from  morpholino  knockdown  of  nfl  a  and  nflb  in  zebrafish 
(Padmanabhan  et  al.,  2009;  Lee  et  al.,  2010).  These  defects  were 
observed  at  48  and  72  hpf  and  probably  resulted  from  impairment 
of  both  maternal  and  zygotic  gene  expression.  In  our  stable 
compound  mutants,  we  observed  partially  penetrant  insufficiency 
of  the  atrioventricular  valve  at  3  dpf  along  with  significant  edema 
and  impaired  blood  circulation  associated  with  irregular  heart  rates 
beginning  at  5-6  dpf  (data  not  shown).  However,  we  could  not 
determine  whether  these  effects  were  primary  or  secondary.  It  is 
possible  that  impaired  cardiac  function  is  the  cause  of  death  of  these 
larvae.  The  absence  of  earlier  cardiovascular  phenotypes  is  most 
probably  due  to  the  activity  of  maternal  transcripts  (Abrams  and 
Mullins,  2009);  confirmation  of  this  interpretation  awaits  the 
creation  of  a  maternal  zygotic  mutation. 

Ubiquitous  and  conditional  Nfl  knockout  mice  have  been 
generated  to  investigate  the  role  of  neurofibromin  in  development 
and  tumorigenesis  (Cichowski  and  Jacks,  2001;  Le  and  Parada, 
2007).  Conditional  loss  of  Nfl  with  p53  deficiency  in  mice  results 
in  the  development  of  grade  III  and  IV  astrocytomas  with  full 
penetrance  (Zhu  et  al.,  2005a),  indicating  that  Nfl  mutations  are 
associated  not  only  with  low-grade  but  also  high-grade 
astrocytoma.  Indeed,  two  independent  studies  have  demonstrated 
that  NF1  mutations  are  found  in  about  15-23%  of  human 
glioblastoma  multiformes  (GBMs)  (Parsons  et  al.,  2008;  The 
Cancer  Genome  Atlas  Research  Network,  2008).  Interestingly,  we 
also  observed  gliomas  in  zebrafish  lacking  nfl  and  p53.  Likewise, 
MPNSTs  are  observed  in  both  mouse  and  zebrafish  models.  Thus, 
the  relative  advantages  of  murine  and  zebrafish  systems  can  be 
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leveraged  for  future  studies  aimed  at  developing  therapeutics  for 
these  lethal  complications  of  NF1. 

In  summary,  we  have  developed  and  characterized  zebrafish  lines 
containing  specific  targeted  mutations  in  nfla  and  nflb.  Compound 
deficiency  of  nfla  and  nflb  results  in  lethality  and  predisposes  to 
tumor  formation.  These  studies  provide  a  powerful  new  tool  for 
analysis  of  neurofibromin  function  and  for  the  development  of 
therapies  for  a  common  human  disorder. 

METHODS 
Zebrafish  lines 

The  «/l«A8,  nflb+1°  and  nflb A55  mutant  alleles  were 

generated  by  application  of  modularly  assembled  ZFNs.  The 
«/i«L1247X  mutant  allele  was  generated  by  TILLING.  Our  nfl 
mutant  alleles  were  crossed  into  various  transgenic  lines,  including 
Tg(gfap:GFP)  (Lam  et  al.,  2009),  Tg (soxlO:GFP)  (Thermes  et  al., 
2002;  Carney  et  al.,  2006)  and  Tg {olig2:GFP)  (Shin  et  al.,  2003),  as 
well  as  the p53e7le7  mutant  line  (Berghmans  et  al.,  2005).  Zebrafish 
were  maintained  under  standard  conditions  as  previously  described 
(Westerfield,  2000).  All  experiments  involving  animals  were 
approved  by  the  Institutional  Animal  Care  and  Use  Committees 
of  Harvard  University  and  the  University  of  Pennsylvania. 

TILLING  with  screening  by  CEL-I  method 

Individual  samples  from  a  preconstructed  ‘live  library’  of  pooled 
genomic  DNA  from  ENU-mutagenized  FI  animals  were  used  as  a 
template  for  PCR  with  the  following  PCR  primer  pairs: 
nfla_outer_F,  5'-TGGCAAATAAATGCTGACAGA-3'  and 
nfla_inner_F,  5'-HEX-TTTTTATATCTCATGTTTAGCTCAC- 
AA-3';  nfla_outer_R,  5'-AAGTCTTAAATGGCCTGAGTGG-3' 
and  nfla_inner_R,  5'-6FAM-AAATGGCCTGAGTGGTAATAAA- 
3'.  A  nested  PCR  was  performed  first  using  the  outer  primer  pair 
with  the  following  PCR  conditions:  94°C  for  2  minutes;  30  cycles 
of  94°C  for  30  seconds,  60°C  for  40  seconds  and  72°C  for  1  minute; 
and  72°C  for  5  minutes.  Amplification  was  then  performed  using 
the  inner  primer  pair  with  the  following  PCR  conditions:  94°C  for 
2  minutes;  25  cycles  of  94°C  for  30  seconds,  60°C  for  40  seconds 
and  72°C  for  1  minute;  and  72°C  for  5  minutes.  PCR  products  were 
denatured,  allowed  to  re-anneal,  subjected  to  CEL-I  digestion  and 
separated  by  acrylamide  electrophoresis  using  a  LI-COR  DNA 
analyzer.  Upon  identification  of  a  genomic  DNA  sample  harboring 
a  mutation  in  the  analyzed  region,  the  individual  animals 
comprising  that  genomic  DNA  pool  were  rescreened  to  identify 
the  appropriate  FI  animal  harboring  the  lesion  of  interest.  This  FI 
animal  was  then  outcrossed  to  wild-type  fish  and  progeny  were 
selected  on  the  basis  of  the  presence  of  the  desired  mutation.  A 
genotyping  strategy  was  developed  to  identify  animals  harboring 
the  nonsense  allele  isolated  by  our  TILLING  strategy  using  the 
following  PCR  primers:  nfla_stop_PstI_F,  5'-CTCTCTTCGA- 
CTCTCGCCATCTGCTGTATCAGCTGC-3'  and  nfla_stop_R, 
5'-GAAGCAGAAGGTCATAATCTTGCTGGCTAGGC-3'.  PCR 
conditions  were  95°C  for  2  minutes;  32  cycles  of  94°C  for  30 
seconds,  62°C  for  30  seconds  and  72°C  for  30  seconds;  and  72°C 
for  5  minutes.  This  generates  a  134  bp  PCR  product.  The  wild- 
type  allele  is  resistant  to  PCR  digestion  by  Pstl,  but  the  mutant 
allele  is  not. 
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Modular  assembly  of  ZFNs 

The  desired  three  finger  zinc  finger  proteins  (ZFPs)  were 
constructed  by  a  splice  overlap  extension  PCR  strategy,  with 
individual  finger  modules  amplified  from  a  archive  of  ZFPs  with 
defined  DNA-binding  specificities  (Zhu  et  al.,  2011).  Three 
individual  fingers  (FI,  F2  and  F3)  were  amplified  using  primers 
specific  to  their  desired  backbone  position,  followed  by  an 
overlapping  PCR  step  to  place  them  together  into  a  single  ZFP 
fragment  which  was  cloned,  sequence  verified  and  subsequently 
subcloned  in  frame  with  a  Fokl  nuclease  variant  in  a  pCS2 
expression  plasmid.  The  primers  for  the  three  backbone  positions 
were  as  follows:  FI  forward,  5'-GCGATGG- 

GTACCCGCCCATATGCTTGC-3'  and  FI  reverse,  5'- 
CACTGGAAGGGCTTCTGGCCTGTGTGAATCCGGATGTG- 
3';  F2  forward,  5'-CATCCGGATTCACACAGGCCAGAA- 
GCCCTTCCAGTGTCGCATCTGC-3'  and  F2  reverse,  5'- 
ATGTCGCATGCAAAAGGCTTCTCGCCTGTGTGGGTGCGG 
ATGTG-3';  F3  forward,  5'-CGAGAAGCCTTTTGCATGCGACA- 
3'  and  F3  reverse,  5'-GCGTAGGATCCACCTGTGTGGATCTT- 
GGTGTG-3'.  The  PCR  conditions  for  amplifying  FI  were  98°C  for 
2  minutes  and  15  cycles  of  94°C  for  30  seconds,  68°C  for  30  seconds 
and  72°C  for  20  seconds.  The  PCR  conditions  for  F2  and  F3  were 
98°C  for  2  minutes  and  20  cycles  of  94°C  for  30  seconds,  57°C  for 
30  seconds  and  72°C  for  30  seconds.  The  three  individual  gel- 
purified  PCR  products  for  positions  FI,  F2  and  F3  were  combined 
(15  ng  of  each)  in  a  PCR  using  Advantage  2  HiFi  polymerase 
(Clontech)  and  subjected  to  the  thermal  cycling  of  94°C  for  2 
minutes  and  five  cycles  of  94°C  for  30  seconds,  55°C  for  30  seconds 
and  72°C  for  30  seconds.  Following  this  cycling  program,  FI  forward 
and  F3  reverse  primers  were  added  to  the  reaction  and  thermal 
cycling  was  continued  as  follows:  25  cycles  of  94°C  for  15  seconds 
and  68°C  for  30  seconds.  This  PCR  product  of  approximately  300 
bp  was  gel-purified,  cloned,  sequence  verified  and  subcloned  in 
frame  with  the  DD/RR  or  EL/KK  Fokl  variants  (Miller  et  al.,  2007; 
Szczepek  et  al.,  2007)  with  and  without  the  3'-UTR  of  nanosl 
(Koprunner  et  al.,  2001). 

ZFN  mRNA  injections  and  genotyping  assays 

We  utilized  protocols  similar  to  those  recently  described  (Zhu  et 
al.,  2011).  pCS2-based  expression  plasmids  containing  our 
constructed  ZFNs  were  linearized  downstream  of  the  SV40 
polyadenylation  signal  and  used  as  templates  for  in  vitro 
transcription  of  ZFN  mRNAs  (Ambion).  One-cell  fertilized 
zebrafish  embryos  were  injected  with  varying  amounts  of  DD/RR 
or  EL/KK  Fokl  variant  mRNAs.  Site-specific  ZFN  function  was 
verified  by  a  genotyping  strategy  wherein  introduction  of  mutagenic 
lesions  at  the  target  site  leads  to  loss  of  a  unique  endogenous 
restriction  site.  This  same  strategy  was  later  used  to  genotype  FI 
animals.  For  nfla  and  nflb,  PCR  was  performed  using  the  following 
primer  pairs  and  PCR  conditions:  nfla  F  genotyping  primer,  5'- 
GGTGTGTATGTAAATGGGCTCAATG-3';  nfla  R  genotyping 
primer,  5'-TACAGTTTCCATAAAACCTGACATTTC-3';  nflb  F 
genotyping  primer,  5'-TGCTACCTGCCGGCAGGCTCAG-3'; 
and  nflb  R  genotyping  primer,  5'-ACCTGTGACCATCA- 
TGTTACTGACA-3'.  PCR  conditions  for  nfla  were  94°C  for  2 
minutes;  35  cycles  of  94°C  for  30  seconds,  54.8°C  for  30  seconds 
and  72°C  for  40  seconds;  and  72°C  for  5  minutes.  PCR  conditions 
for  nflb  were  94°C  for  2  minutes;  35  cycles  of  94°C  for  30  seconds, 
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62°C  for  30  seconds  and  72°C  for  40  seconds;  and  72°C  for  5 
minutes.  For  both  nfla  and  nflb,  a  223-bp  PCR  product  was 
generated.  Mutant  alleles  were  resistant  to  subsequent  digestion 
with  BspCNl  or  Ddel  {nfla)  and  Bsll  or  £coNI  {nflb),  whereas  the 
wild-type  alleles  were  completely  digested.  The  molecular  identities 
of  the  various  mutant  alleles  were  determined  by  cloning  and 
sequencing  the  restriction  enzyme-resistant  PCR  products  from 
individual  embryos  derived  from  outcrosses  of  FI  animals  to  wild- 
type  fish. 

Larval  genotyping 

A  modified  fin  clip  genotyping  strategy  was  utilized  to  identify  larval 
genotypes  prior  to  sacrifice.  Briefly,  2-dpf  larvae  were  anesthetized 
with  tricaine  and  stereomicroscopic  amputation  was  performed  on 
the  caudal  fin  with  an  angled  dissecting  knife  (Fine  Science  Tools, 
10056-12).  Genomic  DNA  from  fin-clipped  tissue  samples  was 
prepared  by  collecting  specimens  in  1.5  pi  of  the  surrounding 
medium  and  dispensing  into  PCR  tubes  containing  7.5  pi  of  60 
mM  NaOH,  which  were  incubated  at  95°C  for  20  minutes  followed 
by  4°C  for  5  minutes  with  the  subsequent  addition  of  1  pi  of  1  M 
Tris-HCl  pH  8.  Genotyping  for  nfla  or  nflb  was  performed  as 
described  with  an  increase  in  PCR  cycle  number  to  40.  Genomic 
DNA  prepared  by  this  strategy  was  sufficient  for  a  single  genotyping 
reaction.  Thus,  parental  crosses  were  selected  to  ensure  all  progeny 
were  homozygous  for  the  non-genotyped  nfl  allele  (e.g.  nfla+,~; 
nflb incross).  This  was  verified  by  genotyping  the  homozygous 
mutant  allele  in  sibling  clutchmates. 

Whole-mount  in  situ  hybridization  and  TUNEL  staining, 
immunohistochemistry  and  BrdU  labeling 

Antisense  RNA  probes  were  generated  for  pip  la  (Park  et  al.,  2002) 
and  mbp  (Lyons  et  al.,  2005)  using  a  digoxigenin  RNA  labeling  kit 
(Roche).  A  previously  published  protocol  (Thisse  and  Thisse,  2008) 
was  followed  with  minor  modifications.  Terminal  deoxynucleotidyl 
transferase  dUTP  nick-end  labeling  (TUNEL)  staining  of  larvae  was 
performed  using  the  ApopTag  Peroxidase  In  Situ  Apoptosis  Detection 
Kit  (EMD  Millipore)  with  minor  modifications.  Larvae  were 
subsequently  blocked  in  2%  blocking  reagent  (Roche)  for  4  hours  at 
room  temperature  and  incubated  overnight  in  anti-digoxigenin-POD 
antibody  (Roche;  1:500)  at  4°C.  TUNEL-positive  cells  were  detected 
using  the  TSA  fluorescein  system  (PerkinElmer).  Apoptotic  cells  were 
quantified  by  counting  TUNEL-positive  cells  between  somites  6  and 
15.  For  immunohistochemical  analysis  of  zebrafish  larvae,  we  used 
mouse  anti-BrdU  (Developmental  Studies  Hybridoma  Bank  G3G4; 
1:1000),  mouse  anti-HuC/D  (Invitrogen  A-21271;  1:100),  mouse  anti- 
pERKl/2  (Sigma;  1:200),  mouse  anti-Zrfl  (Zebrafish  International 
Resource  Center;  1:1000),  rabbit  anti-BLBP  (Millipore  AB9558; 
1:1000),  rabbit  anti-pS6  (Cell  Signaling  2211;  1:200),  and  rabbit  anti- 
SOX10  (Park  et  al.,  2005)  (a  generous  gift  from  Bruce  Appel, 
University  of  Colorado,  Denver,  CO;  1:5000)  as  primary  antibodies. 
For  fluorescent  detection  of  antibody  labeling,  we  used  anti-mouse 
and  anti-rabbit  IgG  antibodies  conjugated  with  Alexa  Fluor  488,  568 
and  647  (Invitrogen;  1:200).  Immunohistochemistry  of  adult  zebrafish 
tumor  samples  was  performed  according  to  a  previously  published 
protocol  (Ligon  et  al.,  2004).  Primary  antibodies  included  anti- 
pERKl/2  (Cell  Signaling  4370;  1:200),  anti-pS6  (Cell  Signaling  2211; 
1:50),  anti-phosphorylated  histone  H3  (Cell  Signaling  9706;  1:100), 
rabbit  anti-SOXlO  (Park  et  al.,  2005)  (1:3000),  mouse  anti-GFAP 
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(Sigma  G3893;  1:10,000),  anti-HuC/D  (Invitrogen  A-21271;  1:200) 
and  anti-synaptophys  (Millipore  MAB5258;  1:1000).  Antibody 
binding  was  detected  using  a  diaminobenzidine-peroxidase 
visualization  system  (EnVision+,  Dako).  Mayers  hematoxylin  was 
used  for  counterstaining.  For  BrdU  labeling,  embryos  were  incubated 
in  BrdU  solution  (10  mM  BrdU  in  2%  DMSO)  for  12  hours.  After 
BrdU  incubation,  embryos  were  fixed  with  4%  paraformaldehyde  and 
embedded  in  1.5%  agarose.  Sections  from  embedded  frozen 
specimens  were  immersed  in  2  M  HC1  for  15  minutes  and  processed 
for  immunohistochemistry.  Paraffin  sectioning  followed  by 
hematoxylin  and  eosin  (H&E)  staining  was  performed  at  the  Dana- 
Farber/Harvard  Cancer  Center  Research  Pathology  Core. 

Behavioral  analysis 

Startle  behavioral  experiments  were  performed  on  5-dpf  larvae  raised 
as  previously  described  (Burgess  and  Granato,  2007b).  Larvae  with 
underinflated  swim  bladders  were  excluded  from  behavioral  testing. 
Acoustic  startle  responses  were  elicited  and  measured  as  previously 
described  (Burgess  and  Granato,  2007a;  Wolman  et  al.,  2011),  such 
that  larvae  could  be  tracked  and  analyzed  individually.  All  startle 
stimuli  were  1000  Hz  waveforms  of  3  milliseconds  duration  at  an 
intensity  of  approximately  150  m/second2.  Stimulus  intensity  was 
calculated  by  measuring  the  displacement  of  the  testing  arena  due 
to  vibration.  To  evaluate  short  latency  C-start  (SLC)  behavior,  images 
were  recorded  30  milliseconds  prior  to  and  90  milliseconds  following 
the  delivery  of  the  3  millisecond  acoustic  stimulus.  To  examine 
acoustic  startle  larval  motor  behaviors,  we  captured  video  recordings 
using  a  MotionPro  high-speed  camera  (Redlake)  at  1000  frames  per 
second  with  512X512  pixel  resolution  using  a  50  mm  macro  lens. 
Behavioral  analyses  were  carried  out  with  the  FLOTE  software 
package  (Burgess  and  Granato,  2007b;  Burgess  and  Granato,  2007a). 
Startle  short-term  habituation  was  performed  and  analyzed  as 
previously  described  (Wolman  et  al.,  2011).  Larvae  were  genotyped 
following  behavioral  testing. 

Western  blotting 

Protein  lysates  were  prepared  from  3-dpf  wild-type,  nfla~'~;  nflb +/+, 
nfla+,+;  nflb~'~  and  nfla~,_ ;  nflb~r~  larvae.  Briefly,  groups  of  10- 
20  larvae  with  identical  genotype  were  anesthetized  with  tricaine, 
deyolked  in  a  solution  of  ice-cold  PBS  with  0.1%  Tween-20  (PBST), 
transferred  to  a  pre-chilled  microcentrifuge  tube  containing  5  pi 
of  lysis  buffer  (150  mM  NaCl,  50  mM  Tris-HCl  pH  8,  5  mM  EDTA, 
1  mM  PMSF,  1%  Igepal  CA  630,  0.5%  sodium  deoxycholate,  0.1% 
SDS  and  IX  Halt  protease  inhibitor  cocktail  (Thermo  Fisher 
Scientific)  per  larva,  sonicated  using  a  Bioruptor  (Diagenode)  and 
cleared  by  centrifugation.  Larval  (100  pig)  and  293T  (25  pig)  protein 
lysates  were  separated  by  gel  electrophoresis,  transferred  to  PVDF 
membranes  and  probed  overnight  at  4°C  with  the  following  primary 
antibodies:  anti-neurofibromin  (Abeam  abl7963;  1:1000),  anti-a- 
tubulin  (Sigma  T6074;  1:1000),  anti-pERI<l/2  (Cell  Signaling  4377; 
1:1000),  and  anti-ERI<l/2  (Cell  Signaling  9102;  1:1000).  Primary 
antibody  binding  was  visualized  on  X-ray  film  using  anti-mouse- 
HRP  (Cell  Signaling  7076;  1:10,000)  or  anti-rabbit-HRP  (Cell 
Signaling  7074;  1:10,000)  secondary  antibodies  along  with 
LumiGLO  (Cell  Signaling)  or  SuperSignal  West  Femto  (Thermo 
Fisher  Scientific)  chemiluminescent  substrates.  Each  Western  blot 
was  performed  in  three  independent  replicates  with  a  representative 
image  of  one  presented. 
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Quantification  of  lateral  line  and  head  melanophores 

Lateral  stripe  melanophores,  including  those  observed  along  the 
horizontal  myoseptum,  were  counted  by  a  blinded  observer  in  live 
progeny  from  nfla+,~;  nflb+l~  incrosses  at  3  and  6  dpf.  Each  discrete 
melanized  region  was  counted  as  a  single  melanophore.  Following 
6-dpf  lateral  stripe  melanophore  quantitation,  larvae  were 
immersed  in  E3  medium  containing  5  mg/ml  epinephrine  (Sigma 
E4375)  for  10  minutes  to  induce  contraction  of  melanosomes 
around  cell  bodies  and  allow  evaluation  of  head  melanophore 
numbers.  Larvae  were  fixed  overnight  in  4%  paraformaldehyde  in 
PBS  at  4°C,  washed  twice  with  PBS  for  5  minutes,  and  melanophores 
anterior  to  somite  one  across  the  crown  of  the  head  (supplementary 
material  Fig.  S14)  were  counted.  Larvae  were  subsequently 
genotyped  for  nfla  and  nflb.  Data  analysis  was  performed  by  one¬ 
way  ANOVA  with  Dunnett’s  post-test  (GraphPad  InStat  3.1a, 
GraphPad  Software). 

Serial  tracking  of  lateral  line  melanophores 

Lateral  stripe  melanophores  of  live  progeny  from  nfla+,~;  nflb or 
wild-type  incrosses  were  individually  imaged  at  3  dpf  followed  by 
incubation  with  0.2  mM  Af-phenylthiourea  (PTU,  Sigma)  to  prevent 
melanin  synthesis.  At  5  dpf,  lateral  stripe  melanophores  were  again 
imaged,  after  which  PTU  was  washed  out.  Larvae  were  reimaged  at 
6  dpf  and  subsequently  genotyped  for  nfla.  Images  were  acquired 
using  a  Nikon  SMZ1500  microscope  and  NIS-Elements  F2.20 
software  with  identical  settings.  Melanophores  in  a  200-pm  region, 
corresponding  to  roughly  12  somites  at  3  dpf,  were  counted  at  3  and 
6  dpf  for  each  larva.  Statistical  analysis  was  performed  using  a  one- 
tail,  unpaired  1-test  (GraphPad  Prism  5,  GraphPad  Software). 

Intestinal  transit  assays 

Groups  of  5-dpf  wild-type  and  nfla~!~;  nflb  larvae  were  placed 
in  individual  wells  of  six- well  plates  containing  feeding  medium  [4 
ml  E3,  2  ml  paramecia  culture  and  5  pi  of  2 -pm  yellow-green 
microspheres  (Polysciences  18338)].  Larvae  were  incubated  at 
28.5°C  for  1  hour  followed  by  several  E3  washes.  Individual  larvae 
were  transferred  to  24-well  plates  and  visually  assessed  for  the 
presence  of  fluorescent  microspheres  in  the  intestinal  bulb.  Extent 
of  intestinal  transit  was  observed  at  2,  4,  6  and  24  hours.  Transit 
was  considered  to  be  complete  when  fluorescence  was  no  longer 
detectable  in  the  intestinal  tract. 

Transmission  electron  microscopy 

TEM  was  carried  out  at  the  Harvard  Medical  School  Electron 
Microscopy  Facility.  Briefly,  embryos  were  fixed  in  2%  formaldehyde 
and  2.5%  glutaraldehyde,  post-fixed  with  1%  osmium  tetroxide  and 
embedded  in  epon.  Sections  were  collected  in  the  trunk  region  of 
embryos.  Images  were  captured  by  a  Tecnai  G2  Spirit  BioTWIN 
electron  microscope  with  an  AMT  2k  CCD  camera. 

Tumor  identification 

«/iaA8/+;  nflb+10l+1°;  p53e7,e7  fish  were  incrossed  and  progeny  were 
manually  evaluated  weekly  for  45  weeks.  Animals  identified  as  having 
tumors  were  separated  and  fin-clipped  for  genotyping  purposes. 
These  samples  were  subsequently  subjected  to  histological  and 
immunohistochemical  analyses  as  described  above  for  determination 
of  tumor  type.  At  the  completion  of  45  weeks,  all  tumor-free  fish 
were  genotyped  for  subsequent  statistical  analysis. 
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Clinical  issue 

Neurofibromatosis  type  1  (NF1)  is  one  of  the  most  commonly  inherited  human 
genetic  disorders.  Despite  nearly  complete  penetrance,  the  clinical  expression 
of  NF1  varies  widely,  even  within  families  harboring  identical  mutations  at  the 
NF1  locus.  Not  surprisingly,  few  genotype-phenotype  relationships  have  thus 
far  been  reported  for  NF1,  suggesting  important  contributions  from  unlinked 
modifier  genes  and/or  environmental  factors.  However,  these  observations  do 
not  preclude  mutations  or  deletions  within  the  NF1  locus  from  influencing 
pathology.  Instead,  they  highlight  the  need  for  better  experimental  tools  to 
address  this  important  and  clinically  relevant  observation.  Additional  models 
of  NF1  are  needed  to  begin  elucidating  these  mechanisms  using  scalable 
chemical  and  genetic  approaches. 

Results 

The  authors  employ  zinc  finger  nucleases  and  TILLING  to  isolate  null  alleles  of 
the  two  zebrafish  orthologues  of  human  NF1,  nfla  and  nflb.  They  report  that 
zebrafish  lacking  nfla  and  nflb  exhibit  valvular  insufficiency,  defects  in 
learning  and  behavior  and  early  larval  lethality.  Larvae  carrying  a  single  wild- 
type  allele  of  either  nfla  or  nflb  are  viable  and  fertile,  suggesting  functional 
redundancy.  The  authors  also  observe  hyperplasia  and  aberrant  differentiation 
in  the  oligodendrocyte  progenitor  cells  and  Schwann  cells  populating  the 
nervous  systems  of  nfl-null  larvae.  This  is  accompanied  by  irregularities  in  the 
myelin  sheaths  surrounding  the  neuronal  axons  of  the  central  nervous  system. 
Human  NF1  is  a  potent  tumor-suppressor  gene  and  the  authors  provide 
evidence  that  zebrafish  nfla  and  nflb  function  similarly:  they  demonstrate 
that  Ras  is  hyperactivated  in  the  spinal  cords  of  nff-null  larvae,  and  that  the 
combined  loss  of  nfl  and  p53  accelerates  tumorigenesis.  Finally,  the  authors 
characterize  a  melanophore  defect  resulting  from  nfl  loss  that  disrupts  the 
uniform  pigmentation  pattern  observed  along  the  lateral  stripes. 

Implications  and  future  directions 

Using  zebrafish  to  probe  the  genetic,  epistatic  and  environmental  factors 
underlying  NF1  pathology  offers  several  important  advantages  over  currently 
available  murine  models.  The  low  costs  and  high  fecundity  of  zebrafish 
coupled  with  their  ability  to  survive  for  several  days  as  haploid  organisms 
make  them  amenable  to  large-scale  genetic  screens.  Thus,  nfl-deficient 
zebrafish  should  greatly  facilitate  the  identification  of  modifier  genes 
influencing  NF1  pathogenesis.  In  addition,  genetic  rescue  experiments  using 
specific  NF1  mutations  or  deletions  could  clarify  the  molecular  basis  of 
pathology.  The  feasibility  of  high-throughput  chemical  screening  using  this 
model  should  provide  additional  valuable  mechanistic  insights  and  identify 
lead  compounds  for  future  therapeutics.  Few  treatment  options  are  currently 
available  for  individuals  affected  with  NF1,  so  advances  in  this  area  are 
urgently  needed.  Importantly,  this  represents  the  first  animal  model  that 
demonstrates  pigmentation  defects  analogous  to  the  pathognomonic  cafe-au- 
lait  spots  seen  in  affected  individuals.  Therefore,  this  model  will  provide  a 
platform  for  further  investigation  of  one  of  the  most  common  clinical 
pathologies  associated  with  NF1. 
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